AD-A035  737 


unclassified 

I OF  3 
AD  A 
035  737 


CORNELL  UNIV  ITHACA  N Y DEPT  OF  OPERATIONS  RESEARCH  F/G  21/5 

SPAERS;  SIMULATION  FOR  THE  PERFORMANCE  AIRCRAFT  ENGINE  REPAIR  S— ETC(U) 
MAY  76  H GIVRAY.  R SLON  N00014-75-C-1172 


U.S.  DEPARTMENT  OF  COMMERCE 
National  Technical  Information  Service 

AD-A035  737 


SPAERS:  SIMULATION  FOR  THE  PERFORMANCE 
OF  AIRCRAFT  ENGINE  REPAIR  SYSTEMS 


Cornell  University^  Ithaca^  New  York 


May  1976 


J 


COLLEGE  OF  ENGINEERING 
CORNELL  UNIVERSITY 

ITHACA,  NEW  YORK  14853 


REPRODUCED  BY  _ . . 

NATIONAL  TECHNICAL 

information  service 

U.S.  DEPARTMENT  OF  COMMERCE 
SPRINGFIELD,  VA.  22161 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  leleosej 
Distribution  Unlimited 


DEPARTMENT 

OF 

OPERATIONS  RESEARCH 


SCHOOL  OK  OPERATIONS  RESEARCH 
AND  INDUSTRIAL  ENGINEERING 
COLLEGE  OF  ENGINEERING 
CORNELL  UNIVERSITY 
ITHACA,  NEW  YORK 


Master  of  Engineering  Report 


May  197-^ 


SIMULATION  FOR  THE  PERFORMANCE 
OF  AIRCRAFT  ENGINE 
REPAIR  SYSTEMS 


by 


Henry  S.  Glvray  and  Robert  A.  Sion 


t... 


I 

I 


S T 


Prepared  for  Naval  Weapons  Engineering  Support  Activity,  Weapons 
Systems  Analysis  Department  (ESA-19),  Washington  Navy  Yard, 
Washington,  D.C.  2037^,  Task  NR  042-335- 


Copy  to  DOC 

Ivlly 


DISTBlBUnON  gfATFMrNT  A~ 


Apptored  fox  public  release; 
Distribution  Unlimited 


Q 


D D C 

f?mf?nnnE 


FEB  16  ign 


tttJ  10 

JlkisE 


UiriLT' 


D 


1 

J 


School  of  Operations  Research 
and  Industrial  Engineering 
Upson  Hall 
Cornell  University 
Ithaca,  New  York  1^853 

May  28,  1976 


Dr.  James  Matthesen 

Naval  Weapons  Engineering  Support  Activity 
Weapons  Systems  Analysis  Department  (ESA-19) 
Washington  Navy  Yard 
Washington,  D.  C.  20374 


Dear  Dr.  Matthesen; 

It  is  our  pleasure  to  submit  this  report  and  documentation  on 
the  Simulation  for  the  Performance  of  Aircraft  Engine  Repair 
S^ystems  (SPAERS)  to  you.  It  was  designed,  and  fully  meets  the 
requirements  for  the  degree  of  Masters  of  Engineering  at  the 
School  of  Operations  Research  and  Industrial  Engineering  at 
Cornell  University. 

We  trust  that  it  will  be  of  service  to  you  and  your  staff  in 
verifying  and  comparing  Naval  aircraft  engine  repair  system 
configurations . 


Sincerely, 


Henry  S.  Givray 


Robert  A.  Sion 


ai'jcx^ 


/ I 
oijOu 


C(H\^ 


OOVt  AT  Cl 


<ii  iitMoRi  a fi  111 

;-5earch  Hepor 


'3VI\liH3:  SIMULATION  FOE?  TIEE^  PLEtl-’ORMANCE  OF 

AIRCRAFT  ENGINE  REPAIR  SYSTEMS 


/.  Al  T nO^(9) 

Henry  Givray 
Robert  Sion 


N0001i4-75-C-1172 


JCRrOPMINCOHOAMiZATlON  NAME  AND  ADORE  SS 

ScEiool  of  Operations  ResearcEi  and 
Ind.  Eng. 

Cornell  Univ.  , Itliaca,  NY  14853 


12.  REPORT  CATE 

May  1976  ■ 


. CCri  T ROLLING  OPFiCE  NAME  AND  ADDRESS 

Office  of  Naval  ResearcEi  (code  436) 
Dept,  of  the  Navy 

Arlington,  VA  22217  

;;  CORTROEUING  OFFICE  NAME  AND  ADDRESS 


L * SSI  FICATICN/  DOWN  Ding 

ECULE 


01S7  Ri  5JT|0N  ST  AT  £mCN  7 (of  (hi a Feport) 


distribution  unlimited 


Approved  for  public  releas 


UISTR'SUTjON  STA7cMENT  (ef  th  9 mbftrect  rnfttmd  In  Stock  20,  II  /root  Report} 


. SURRLEMEkTARY  NOTES 

Also  supported  by  the  Naval  Weapons  Engineering  Support  Activity 
(ESA-19)  - Systems  Analysis,  Washington  Navy  Yard,  Washington,  D 


'Cnatcjy  •artfl  Itien’.liy  by  block  nu:rb9i) 


KEY  WORDS  fConnnoo  on  rrver^•  aid 

Aircraft  Engine 
Simulation 

Multi-echelon  System 
Navy  Supply  System 


The  allocation  of  spare  aircraft  engines  is  critical  to  the 
U.S.  Naval  aircraft  operation’s  performance.  An  aircraft  in  this 
system  becomes  Inoperative  in  the  event  of  an  engine  failure  and 
remains  in  that  state  until  it  is  replaced  by  a serviceable  engine 
An  engine  is  removed  upon  failure  and  subsequently  is  recovered 
by  repairing  it  at  the  location's  repair  facilities  or  elsewhere. 
However,  the  availability  of  a spare  engine  at  the  location  could 
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r’<-fiuci'  tho  amcunt  of  timo  whichi  an  aircraft  zocnt  in  an  Inoperative 
:;tato  due  to  (•rrf’;Lne  rcfiair  timo. 

Anal  y t j c:i  i model:;  have  been  dovol(;pod  by  the  DOD  to  calculate' 
c.pan.'  enp. Ine  r-equ  i x-eiMorit.a  throu/^hout  the  ayr-toin.  The  model:;  roqulrr 
.aa.-.nmpt  i on;;  Ic  'iiade  about  (;c;rtaln  partr;  c;f  the  ayatem.  One  lmpo:-tan 
a ,;;-.ii(i,f)t  Ion  i :;  that  carrier:;  ar<;>  rnaaimed  to  be  continuously  connect<rJ 
1.0  the  ropai./'  ;;,ystoin.  ThU;  a:;t.umpt  i on  Iprior-.-s  the:  fact  that  in 
[•f  ility  tii'M.r  only  connection  to  the  repair  system  Is  at  a port.  Tii 
i ;ri[)l  Ication  is  ll;at  under  this  ac.r.umption,  carx-ic't'S  may  I'.v.aap  a fa.lle 
en.",jne  for  a serviceable  one  at  any  time,  dii'octly  wltli  .th.e  repair 
1 ucat Ion . 

In  Lite  re.al  situations,  carriers  perfor.m  this  trade  usinr,  the 
p-oi't  as  an  intermediary  and  only  at  Its  ti.me  of  Its  docking.  Trie 
effect  is  that  in  reality,  ports  see  mor-e  engines  in  their  stock  tii-a 
they  do  in  analytical  models.  If  the  port  supports  a flying 
a-'t-ivlty  at  that  location,  the  additional  engine  availability  increa 
it.s  flying  performance  and  may  decrease  that  of  the  carriers. 

A Simulation  for  the  Performance  of  Aircraft  Engine  Repair 
System  (SPAERS)  was  developed  to  simulate  different  configurations 
of  an  aircraft  repair  system.  The  analysis  section  of  this  report 
shows  a comparison  between  the  two  situations  described  above,  namel 
the  analytical  rendering  of  tiie  r’epalr  system,  and  a situation  more 
closely  resembling  the  real  syste.m  dynamics. 

To  a significant  level  (.05),  differences  in  the  performance 
criterion*  were  detected  between  the  two  configurations  of  the  same 
system  at  the  port  whch  supported  flying  activity,  and  on  one  of  the 
more  active  carriers. 

Further  study  of  similar  repair  systems  should  be  made  to  suppo 
inoro  general  statements  and  observations  concerning  the  relations 
of  carriers  and  ports  to  the  entire  system. 

It  was  shown,  however,  that  SPAERS  can  be  used  as  a tool  for 
verifying  the  performance  predicted  by  analytical  models  for 
:;eiocted  repair  systems,  and  enables  thier  subsefpi;nt  statistical 
co:!ipar i son . • 


(.'riterion  is  defined  as:  average  down  aircraft  at  a base  divided 

by  number  of  aircraft  at  that  base.  A down  aircraft  is  defined  as 
an  aircraft  unable  to  'fly'  due  to  the  failure  of  one  of  its  engines 
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PREFACE 


In  any  system  where  the  availability  of  physical 
components  determines  the  operational  quality  of  the  system, 
inventory  management  plays  a key  role.  The  United  States 
Navy's  aircraft  operation  is  an  example  of  this  situation. 

The  physical  components  in  this  case  are  aircraft  engines, 
and  the  operational  quality  of  the  system  is  measured  by 
aircraft  readiness.  Here  the  components  fail  over  time  and 
can  be  subsequen'^ly  recovered  through  the  repair  system.  A 
failed  engine  renders  an  aircraft  inoperative  until  a 
serviceable  engine  becomes  available  as  a replacement.  When 
a demand  for  an  engine  is  not  met,  it  is  said  to  be  back- 
ordered and  met  at  a later  date.  To  compensate  for  the  time 
associated  with  repairing  an  engine,  spare  engines  are 
pre-positioned  at  various  locations.  These  locations  which 
support  repair,  supply,  and  flying  activity  are  called  bases. 

The  Department  of  Defense  presently  uses  a method, 
described  in  DODI  ^<230.4  [1],  to  establish  spare  engine 
stock  levels  for  each  base.  The  method's  objective  is  to 
stock  each  base  so  that  the  base  has  outstanding  backorders 
for  only  a small  proportion  of  time. 

The  DODI  ^<230. 4 method  of  establishing  spare  stock 
levels  at  a particular  base  does  not  consider  an  Important 
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system  Interaction.  In  particular,  it  does  not  account  for 
the  manner  in  which  resupply  times  may  vary  with  the  spare 
stock  levels  at  resupply  bases. 

Consequently,  a model  called  NAVMET  [2]  was  developed 
which  explicitly  considers  the  effect  of  stocking  inventory 
at  one  location  on  the  availability  of  engines  at  other  loca- 
tions. NAVMET' s objective  is  to  optimize  the  quantity  and 
allocation  of  spare  engine  with  respect  to  a constraint  on 
average  down  aircraft.  It  is  applied  in  the  environment  of 
the  Navy's  multi-echelon  repair  system. 

The  dynamic  nature  of  the  Navy's  repair  system, 
particularly  the  aircraft  carriers'  relation  to  the  system 
have  raised  some  questions  concerning  both  models'  assumptions 
about  the  demand  process.  These  questions,  which  will  be 
developed  more  fully  in  tnis  report,  formed  the  basis  for  a 
Masters  of  Engineering  at  Cornell  University.  The  project 
involves  modelling  the  dynamics  of  the  actual  system  by 
means  of  a digital  computer  simulation.  The  simulation  will 
be  used  to  study  the  effect  of  an  assumption  made  in  the 
analytic  models  regarding  carrier  operations  on  expected 
system  performance. 

This  simulation,  called  SPAERS  (Simulation  for  the 
Performance  of  Aircraft  Engine  Repair  S^ystems),  was  developed 
for  the  Naval  V/eapons  Engineering  Support  Activity  (ESA-19). 

It  was  designed  as  a tool  for  observing  the  long  run 
performance  of  a repair  system  given  its  system  parameters. 
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SPAERS  input  may  be  adjusted  to  simulate  the 
environment  in  either  the  configuration  assumed  by  the 
analytical  models,  or  the  one  which  more  closely  resembles 
the  real  system.  This  will  enable  a comparison  of  the 
system's  performance  to  be  made  between  the  two  modes  of 
operation,  given  the  same  spare  engine  stock  levels. 

The  authors  wish  to  thank  the  faculty  and  members  of 
the  School  of  Operations  Research  and  Industrial  Engineering 
at  Cornell  University  for  their  valuable  consultation.  In 
particular:  Professors  John  Muckstadt,  William  Maxwell  and 

Thomas  Santner;  and  Dr.  James  Matthesen  and  Mr.  William 
Booth  of  the  Naval  Weapons  Engineering  Support  Activity. 
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ABSTRACT 

The  allocation  of  spare  aircraft  engines  is  critical 
to  the  U.S.  Naval  aircraft  operation's  performance.  An 
aircraft  in  this  system  becomes  inoperative  in  the  event  of 
an  engine  failure  and  remains  in  that  state  until  it  is 
replaced  by  a serviceable  engine.  \n  engine  is  removed  upon 
failure  and  subsequently  is  recovered  by  repairing  it  at  the 
location's  repair  facilities  or  elsewhere.  However,  the 
availability  of  a spare  engine  at  the  location  could  reduce 
the  amount  of  time  which  an  aircraft  spent  in  an  inoperative 
state  due  to  engine  repair  time. 

Analytical  models  have  been  developed  by  the  DOD  to 
calculate  spare  engine  requirements  throughout  the  system. 
The  models  require  assumptions  be  made  about  certain  parts 
of  the  system.  One  Important  assumption  is  that  carriers 
are  assumed  to  be  continuously  connected  to  the  repair 
system.  This  assumption  Ignores  the  fact  that  in  reality 
their  only  connection  to  the  repair  system  is  at  a port. 

The  implication  is  that  under  this  assumption,  carriers  may 
swap  a failed  engine  for  a serviceable  one  at  any  time, 
directly  with  the  repair  location. 

In  the  real  situations,  carriers  perform  this  trade 
using  the  port  as  an  intermediary  and  only  at  its  time  of 
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its  docking.  The  effect  is  that  in  reality,  ports  see  more 
engines  in  their  stock  than  they  do  in  analytical  models. 

If  the  port  supports  a flying  activity  at  that  location,  the 
additional  engine  availability  increases  its  flying  per- 
formance and  may  decrease  that  of  the  carriers. 

A Simulation  for  the  P^erformance  of  Aircraft  Engine 
R^epair  S^ystem  (SPAERS)  was  developed  to  simulate  different 
configurations  of  an  aircraft  repair  system.  The  analysis 
section  of  this  report  shows  a comparison  between  the  two 
situations  described  above,  namely  the  analytical  rendering 
of  the  repair  system,  and  a situation  more  closely  resembling 
the  repl  system  dynamics. 

To  a significant  level  (.05),  differences  in  the 
performance  criterion*  were  detected  between  the  two 
configurations  of  the  same  system  at  the  port  which  supported 
flying  activity,  and  on  one  of  the  more  active  carriers. 

Further  study  of  similar  repair  systems  should  be 
made  to  support  more  general  statements  and  observations 
concerning  the  relations  of  carriers  and  ports  to  the  entire 
system . 

It  was  shown,  however,  that  SPAERS  can  be  used  as  a 
tool  for  verifying  the  performance  predicted  by  analytical 


*Criterlon  is  defined  as:  average  down  aircraft  at 

a base  divided  by  number  of  a.ircraft  at  that  base.  A down 
aircraft  is  defined  as  an  aircraft  unable  to  ’fly'  due  to 
the  failure  of  one  of  its  engines. 
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INTRODUCTION 

United  States  Naval  aircraft  operations  depend  on 
the  quantity  and  distribution  of  spare  aircraft  engines 
throughout  Its  aircraft  engine  repair  systems.  As  aircrafu 
fall  over  time,  the  failed  engines  are  subsequently  recovered 
by  the  repair  system.  Aircraft  down  time  due  to  repair 
system  lead  times  can  be  I’educed  by  allocating  spare  engines 
to  particular  system  locations  (bases).  An  engine  failure 
is  considered  a demand  upon  the  system  which  can  only  be 
satisfied  by  a serviceable  engine.  If  a demand  cannot  be 
satisfied  immediately,  it  is  backordered  and  filled  as  soon 
as  possible. 

One  method  of  spare  engine  allocation  presently  used 
by  the  Navy  is  described  in  the  Department  of  Defense 
Instruction  4230.^1  [1].  Its  objective  is  to  stock  each  base 
for  a small  proportion  of  time.  This  method,  however, 
ignores  the  effect  of  the  spare  engine  supply  levels  on 
resupply  times. 

NAVMET  [2]  takes  stock  dependent  resupply  times  into 
consideration . Its  objective  is  to  maintain  the  expected 
num.ber  of  backorders  to  flying  activity  at  specified  minimal 
levels.  MAWET  prouuces  an  optimal  solution  to  the  problem 
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of  allocating  quantities  of  spare  engines  within  the  Navy’s 
four  echelon  repair  systems  with  respect  to  a constraint  on 
average  time  weighted  downed  aircraft. 

The  system  dynamics,  particularly  those  Involving 
aircraft  carriers,  present  some  reasonable  doubt  as  to  the 
validity  of  the  demand  process  that  is  assumed  for  carriers 
in  both  of  these  models. 

The  following  report  presents  the  development  of  a 
computer  simulation  model,  called  SPAERS  (Simulation  for  the 
Performance  of  Aircraft  Engine  Repair  S^ystems),  which  serves 
to  statistically  compare  the  actual  operation  of  the  repair 
system  with  the  manner  in  which  it  is  assumed  to  operate  by 
the  analytical  models.  The  purpose  of  this  report  is  to 
present  the  system  characteristics  which  are  relevant  to 
model  formulation,  a discussion  of  the  analytical  methods 
for  establishing  stock  levels  in  the  system,  and  the  moti- 
vation for  a computer  simulation  of  the  system.  As  SPAERS 
will  be  used  by  the  Naval  V/eapons  Engineering  Support 
Activity,  as  a selective  verification  procedure  for  predicted 
system  performance,  the  manner  in  which  the  simulation 
treats  the  real  system  is  discussed  in  detail.  This  dis- 
cussion includes  the  simulation's  initial  conditions,  model 
assumptions,  and  user  availabilities. 

Analysis  Includes  the  results  of  a specific  system 
simulation.  Stock  levels  were  computed  for  the  system  using 
the  DODI  calculation,  and  two  runs  were  performed.  One 
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operated  the  system  In  the  manner  which  is  assumed  by  the 
I model,  and  the  second  run  operated  the  system  using  the 

actual  system  dynamics.  Finally,  some  observations  were  made 
concerning  the  role  of  the  port  in  computation  of  stock 
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RELEVANT  SYSTEM  CHARACTERISTICS 

In  order  to  aid  the  discussion  of  the  repair  system 
in  later  sections  of  this  report,  some  features  of  the 
aircraft  engine  repair  environment  which  are  important  and 
often  assumed  by  a mathematical  formulation  will  be  described. 

A base  is  defined  as  a location  where  supply,  repair, 
or  flying  activity  occurs.  The  base’s  flying  activity  gener- 
ates failed  engines  which  are  sent  to  the  repair  function. 

If  the  repair  function  at  this  base  does  not  have  the 
facilities  to  perform  the  repair  required,  the  engine  is 
sent  to  another  repair  base  in  exchange  for  a serviceable 
(operational)  engine. 

Engines  which  complete  repair  at  a base  become  part 
of  the  base’s  serviceable  stock  supply.  The  supply  function 
is  responsible  for  maintaining  the  flying  activity  at  its 
own  base.  In  addition,  it  must  exchange  a serviceable 
engine  for  a failure  which  is  sent  to  its  own  location  for 
repair  (see  figure  1). 

The  base  repair  facility’s  ability  to  deal  with 
failures  depends  on  the  sophistication  of  its  repair 
equipment.  A base  is  identified  by  the  most  complex  repair 
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BASE  ACTIVITY 
(engine  flows) 


that  it  can  accomplish. 


Five  degrees  describe  engine  repair 


no  repair 
3®  minor  repair 
2°  moderately  complex  repair 
1®  major  repair 
0®  overhaul*;  major  repair 

If  an  engine  falls  at  a base  which  cannot  perform 
the  required  degree  of  repair,  the  engine  is  sent  to  a 
higher  echelon  (lower  degree)  repair  base.  In  this  case  the 
flying  activity  at  which  the  failure  occurred  demands  an 
engine  from  the  supply  function  at  its  base.  The  fact  that 
an  engine  was  sent  for  repair  to  another  base  requires  that 
the  repair  base  return  a serviceable  engine  to  the  supply  at 
the  base  which  experienced  the  failure  (see  figure  2).  The 
repair  system  is  fixed  in  the  sense  that  all  J®  failures 
which  occur  at  base  I are  always  repaired  at  a base  K for 
all  I and  ^*». 

A supply  function's  inability  to  meet  any  demands  is 
registered  as  a backorder  and  filled  as  soon  as  possible. 

A repair  system  (see  figure  3)  is  defined  as  a set 
of  bases  having  a common  0®  repair  base  for  a particular 
engine  model  type.  The  types  of  bases  which  comprise  a 
system  are  described  as  follows. 


•Periodically,  engines  are  removed  from  aircraft  for 
scheduled  overhaul  even  though  they  have  not  failed. 

••Occasionally , in  the  real  system,  a 2®  failure 
occurring  at  a 2®  repair  base  will  be  sent  to  a 1°  repair 
base.  A similar  situation  occurs  with  1°  failures  and 
overhaul  locations. 
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DEHANDS 


A (2®) 


demand 

R repair 
FA  flying  activity 
S supply 


In  the  above  figure,  base  A is  a 2°  repair 
base  and  B Is  a 1°  base.  Suppose  the  flying  activity 
at  A has  experienced  a 1°  failure.  Thus  A's  repair 
facility  cannot  repair  the  engine,  and  It  must  be 
sent  to  B.  The  two  demands  which  occur  are: 

1)  Flying  activity  on  Supply  Function  at  A,  and 

2)  Supply  Function  at  A on  Supply  Function  at  B. 
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1)  A carrier  has  Its  own  "on  board"  supply  function 
and  flying  activity.  When  a carrier  is  deployed,  its  link 
with  higher  echelon  repair  is  broken  and  its  flying  activity 
performs.  A non-deployed  carrier  is  connected  to  the  repair 
system,  but  its  flying  is  inactive. 

Due  to  the  carrier’s  particular  relationship  with 
the  rest  of  the  repair  system,  on-board  engine  failures 
which  the  carrier  is  unable  to  repair  will  not  enter  the 
repair  system  until  the  carrier  becomes  non-deployed 
(docked)*  . 

During  the  dock  time,  the  port  sends  the  carrier's 
failed  engines  to  their  respective  repair  bases  in  exchange 
for  serviceable  engines.  If  some  carrier  demands  are  not 
satisfied  durlngs  its  docking,  they  are  backordered  at  the 
port . 

2)  A port  has  full  responsibility  for  supplying 
carriers  before  they  are  deployed.  A port  may  have  any 
number  of  carriers  to  service,  but  a carrier  is  always 
serviced  by  one  port.  The  port  may  also  support  its  own 
flying  activity. 

3)  The  depot  performs  overhauls  (0°  repair)  for  the 
system  and  does  not  generally  support  its  own  flying  activity. 

The  "pipeline"  describes  the  connections  between 
bases  or  intra-base  functions  and  is  measured  in  time 
durations.  Engines  may  be  said  to  flow  through  the  pipeline 


1 


*In  special  cases  in  the  acutal  system,  a serviceable 
engine  may  be  delivered  to  a carrier  while  it  is  deployed. 
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at  a constant  rate. 

A "twisted  pipeline"  is  used  to  explain  the 
connection  between  carriers  and  port.  When  the  carrier  is 
deployed,  its  pipeline  to  the  port  is  cut  off,  or  in  a sense, 
"twisted"  so  that  the  flow  of  engines  may  not  reach  the  port 
until  the  carrier  docks.  It  may  be  noted  here  that  present 
models  for  determining  spare  engine  requirements  at  a base 
ignore  the  twisted  pipeline  configuration.  Rather,  they 
assume  a fixed,  continuous  resupply  time  for  a carrier  (as 
if  it  were  a ground  base).  The  term,  "continuous  pipeline" 
will  be  used  to  designate  the  treatment  of  carrier  activity 
assumed  by  the  analytical  models.  (The  motivation  for  this 
study  is  the  validation  of  the  continuous  pipeline 
assumptions . ) 

Resupply  time  is  defined  to  be  the  average  amount  of 
time  v/hlch  it  takes  for  a failed  engine  to  be  replaced  by  a 
serviceable  one  at  each  base. 

A base  is  said  to  have  custody  of  an  engine,  in  the 
sense  of  ownership,  if  the  engine  is  en  route  to  the  base, 
or  in  the  base’s  repair  or  supply  activity.  Only  serviceable 
engines  may  be  used  to  satisfy  a demand. 


Ill 


SYSTEM  PARAMETERS 

The  previous  section  described  system  activity  which 
is  relevant  to  mathematical  formulation.  This  section  will 
explain  the  parameters  which  describe  the  repair  system 
operation . 

One  engine  model  type  is  a necessary  constraint  in 
defining  a system.  Relevant  flying  activity,  then,  consists 
of'  those  aircraft  types  having  the  system  engine.  These 
aircraft  are  considered  a base's  permanent  property. 

The  combination  code  describes  the  aircraft  type's 
use  of  the  system  engine.  It  is  associated  with  failure 
rates  and  the  probability  of  J°  failure  for  the  system  engine 
where  J=0,  1,  2,  3. 

The  various  rates,  probabilities  and  times  which  are 
used  as  parameters  for  the  system  can  be  found  in  computer 
output  from  the  U.  S.  Navy's  PASER* . The  system  "geography" 
is  described  by  a matrix  N where  n^^  is  the  base  number  to 
which  a ^ degree  failure  is  sent  given  that  it  occurs  at 
base  N is  easily  derived  from  the  PASER  data. 


- *ProJectlon  of  Aircraft  Spare  Engine  Requirements 

I is  a set  of  computer  software  used  by  the  U.  S.  Navy  for 

‘ computing  Spare  Engine  Requirements. 

! 
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A specific  flying  rate  (in  hours  per  month)  is 
assigned  to  aircraft  types  at  a base.  Any  number  of  aircraft 
types  and  numbers  of  that  aircraft  type  may  be  found  at  a 
base.  A dally  flying  rate  is  assumed  to  be  1/30  of  the 
flying  nours  per  month. 

Travel  time  from  base  ^ to  base  J is  assumed  to  be 
constant  as  are  repair  times  for  a given  failure  degree 
regardless  of  where  that  repair  takes  place. 

The  discrete  probability  distributions  describing 
the  probability  of  a J°  failure  are  available  from  the  same 
PASER  computer  records  as  are  mean  times  between  removal  for 
a given  combination  code. 

Other  system  Inputs  Include  the  number  of  engines 
per  aircraft  (UPA)  and  the  desired  performance  criterion. 

The  criterion  will  be  deflnea  by  the  average  time  weighted 
downed  aircraft  which  are  due  to  engine  failures  divided  by 
the  total  number  of  aircraft  which  utilize  the  system  engine 
at  that  base. 
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IV 

BACKGROUND 

The  system  described  in  the  previous  section  operates 
for  the  purpose  of  maintaining  flying  activity.  Measures 
which  indicate  the  system's  effectiveness  are  in  number  of 
down  aircraft  at  any  point  in  time  and  in  the  durations  of 
inoperative  aircraft.  Other  measures  of  effectiveness  would 
be  described  by  a supply  function's  ability  to  meet  demand, 
and  the  numbers  and  durations  of  its  backorders. 

An  obvious  limitation  on  the  performance  of  a supply 
system  is  the  number  of  spare  engines  that  are  in  the  system 
to  compensate  for  repair  and  travel  times  which  are  Involved 
in  recovering  engines.  At  the  same  time,  it  is  equally 
Important  that  spare  engines  are  located  at  the  "correct" 
supply  function  in  order  to  achieve  acceptable  levels  of 
flying  performance. 

This  section  will  describe  some  of  the  performance 
measures  and  explain  two  models  that  are  used  by  the 
Department  of  Defense  to  allocate  the  appropriate  number  of 
spare  engines  in  an  attempt  to  optimize  using  these  criterion 
measures.  Furthermore,  the  need  for  a simulation  model  to 
statistically  verify  some  of  the  assumptions  made  by  the  DOD 
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models  will  be  demonstrated. 


Measures  of  Performance 

The  ready  rate  may  be  defined  by  the  probability 
that  there  are  less  than  S units  in  resupply  at  a base  where 
S is  a decision  variable. 

S 

RR  = ^p(x|XT), 
x*0 

where  p(x|XT)  is  the  probability  of  having  ^ units  in 
resupply . 

The  major  assumption  here  is  that  demands  follow  a 
Poisson  distribution  with  rate  XT,  where  X = dally  failure 
rate  and  T = average  resupply  time.  Resupply  at  a base  is 
defined  to  be  physical,  on-hand  engines  plus  the  number  of 
engines  in  base  repair  plus  serviceable  stock  from  other 
bases,  less  backorders. 

Fill  rate  is  the  proportion  of  demands  that  can  be 
filled  immediately,  or  the  probability  that  S - 1 units  are 
in  resupply 

S-1 

FR  = J p'xJXT) 
x = 0 

Again  the  demand  process  is  assumed  to  be  Poisson. 

Although  these  measures  are  used  extensively  by  the 
DOD  to  compute  requirements,  they  do  not  reflect  system 
performance  in  a complete  sense.  Neither  the  number  nor  the 
duration  of  backorders  can  be  inferred  from  these  measurements. 

A more  suitable  measurement  for  both  planning  and 
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evaluation  purposes  is  a base's  time-weighted  average 
ji  backorders,  measured  in  backorder-days  per  day.  A mathe- 

i a 

matical  rendering  of  this  measurement  is  given  by: 

; J 00 

B(S)  = I (x-S)p(x|XT) 
x=S+l 

where  S is  the  on-hand  stock  at  the  repair  base.  The 
measure  has  relevance  only  at  bases  with  flying  activity. 
Therefore,  a non-flying  base  would  not  calculate  a backorder 
criterion.  Its  stock,  however,  effects  the  performance  of 
the  flying  bases  it  resupplies. 

This  criterion  can  be  normalized  over  all  bases  by 
dividing  it  by  the  total  aircraft  at  the  base.  Note  that  a 
convenient  aspect  of  this  measurement  is  its  convexity  in  S, 
; given  T,  making  it  a reasonable  objective  in  optimization 

. models . 


Spare  Stock  Level  Determination 

I . 

Given  the  above  performance  measures  and  the  Poisson 
demand  assumption,  the  Department  of  Defense  rationale  for 
supplying  safety  stock  can  be  illustrated.  First,  the 

I 

^ parameters  X and  T are  approximated.  Daily  Demand  Rate  (X) 

for  engines  at  a base  can  be  approximated  by  multiplying 
average  flyji’g  hours  per  month  by  the  number  of  aircraft  at 
the  base,  by  the  failure  rate  per  engine  times  the  average 
number  of  engines  per  aircraft  at  the  base.  This  quantity 
is  in  terms  of  engine  failures  per  day.  The  average  resupply 

I 


time  T is  derived  by  multiplying  the  probability  of  a J° 
failure  at  a base  by  the  travel  time  to  the  J°  repair  base 
and  summing  over  all  ^ plus  the  repair  time  at  base  A times 
the  probability  of  repair  at  base  A (see  figure  ^).  By 
using  the  desired  ready  rate  criterion  of  a,  the  safety 
stock  level  S can  be  calculated  for  the  base  A by  solving 

Sa 

x=0  ^ " 

This,  in  fact,  is  the  method  used  by  the  Department  of 
Defense . 


It  becomes  clear,  that  for  any  value  of  a,  S is 
fixed,  and  the  other  performance  measures  are  implicitly 
determined . 

The  difficulty  with  determining  an  S (spare  stock 
level)  using  the  metnod  above  lies  in  the  fact  that  the 
average  resupply  time  T is  not  truly  represented.  Clearly, 
if  the  repair  base  had  some  on-hand  stock,  resupply  time 
from  a J°  repair  base  could  be  shortened  by  the  time  xt 
takes  to  repair  an  engine.  That  is.  Instead  of  incurring  a 
delay  due  to  repair  time,  an  engine  could  be  sent  off 
immediately.  In  a probabilistic  sense,  this  delay  would 
vary  with  the  amount  of  spare  stock  which  the  repair  base  had. 
T^  might  now  be  represented  as 


+ H(Sj)) 


where 


is  the  amount  of  safety  stock  at  base 


H(Sj)  can 
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be  defined  as  the  number  of  backorder  days  per  demand  such 
that , 

00 

H(S.)  = ^ I (x  - S.)p(x|XT.) 

J x=S+l  ^ J 

where  X = J pjXj  or  the  system  demand  on  the  repair  base. 

Thus,  it  can  be  seen  that  system  interdependence  is 
considered  by  describing  the  resupply  time  as  a function  of 
the  repair  base  stock.  The  NAVMET  [2]  model  employs  the 
above  function  and  develops  an  algorithm  for  determining  the 
, optimal  quantity  and  distribution  of  spare  aircraft  engines, 

i.  The  application  of  this  model  to  a four  echelon 

system  Involves  a dynamic  approach  whereby  subproblems  which 
segregate  bases  into  "families"  are  solved.  Full 
exposition  of  the  algorithm  is  beyond  the  scope  of  this 
[ report.  The  interested  reader  is  referred  to  Naval 

Publication  R-7511  (NAVMET)[2]. 

I The  optimality  of  a solution  obtained  from  the  NAVMET 

;i  algorithm  requires  the  validation  of  some  assumptions  about 

■ i 

the  system.  One  which  requires  discussion  at  this  point, 

f . 

I;  and  serves  ar  the  motivation  for  this  Investigation  concerns 

f ; 

'■  the  continuous  pipeline. 

Recall  that  in  the  actual  system,  aircraft  carriers 
leave  port,  and  for  all  practical  purposes,  are  not  in 
contact  with  the  repair  system  for  the  duration  of  their 
deployment.  In  this  way,  the  repair  pipeline  is  paid  to  be 
twisted . 

1 The  NAVMET  and  the  DODI  ii230.^  approach,  both  assume 
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that  carriers  resemble  ground  bases  and  are  always  connected 
to  the  repair  system  with  a fixed  resupply  time.  Carrier 
resupply  times  are  approximated  by  the  average  resupply  time 
which  they  would  experience  in  the  actual,  twisted  pipeline 
situation . 

Clearly,  the  surges  and  lags  which  the  twisted 
pipeline  carriers  impose  upon  the  repair  system  cast 
aspersious  on  the  ability  of  the  Poisson  process  to  describe 
the  demand  process.  Since  the  description  of  the  demand 
process  is  crucial  to  finding  an  optimal  solution  to  the 
allocation  problem,  it  is  desirable  to  note  the  conditions 
under  which  the  continuous  pipeline  is  a valid  approximation 
to  the  twisted  pipeline. 

The  question,  then,  is  one  of  degree.  On  the  basis 
of  performance  measures,  especially  the  backorder  criterion, 
how  closely  does  the  continuous  pipeline  configuration 
resemble  the  twisted  pipeline?  Other  questions  might  be 
posed;  What  parts  of  a repair  system  feel  the  effect  of  a 
twisted  pipeline?  If  any,  what  relationships  do  carriers' 
continuous  pipeline  performance  have  with  its  demand 
rate  and  deployment  schedules  in  the  twisted  pipeline? 

The  dynamics  of  a twisted  pipeline  renders  it  too 
complex  to  model  analytically.  A digital  computer  simulation 
seemed  to  provide  the  most  convenient  approach  to  an  effec- 
tive comparison.  SPAERS  or  the  Simulation  for  the 
Performance  of  Aircraft  E^nglne  Repair  S^ystems  was  developed 
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to  statistically  compare  performance  of  a twisted  pipeline 
configuration  with  that  of  the  continuous.  In  addition,  the 
NAVMET  performance  predictions  could  also  be  compared  with 
those  of  SPAERS. 


I 


I 


J 


V 


SPAERS 

In  the  previous  section  the  present  methodology  for 
determining  spare  aircraft  engines  was  discussed.  It  was 
pointed  out  that  a procedure  that  would  test  the  validity  of 
the  analytical  models  was  needed.  Furthermore,  It  was 
determined  that  a simulation  approach  could  adequately 
accomplish  this  task. 

The  SPAERS  model  Is  designed  to  examine  the 
performance  of  the  "optimal"  spare  stock  levels  determined 
by  the  anaylltlcal  model,  NAVMET,  under  dynamic  conditions. 

In  this  section  the  SPAERS  model  Is  discussed.  First 
the  underlying  assumptions  and  Initial  conditions  of  the 
model  are  outlined.  SPAERS  options  and  capabilities  along 
with  SPAERS  data  output  are  briefly  described  next.  The 
model's  breakdown  of  the  real-world  system  In  terms  of 
entities  (objects  of  the  system),  their  attributes  (charac- 
teristics of  the  entitles)  and  events  (changes  in  the  state 
of  a system  entity)  is  provided  at  the  end  of  the  section. 

A detailed  description  of  the  construction  and 
functioning  characteristics  of  the  model  may  be  found  in 
Appendix  A.  User  instructions  for  operation  are  provided  in 
Appendix  B. 


21 


22 


A.  Assumptions 


Some  assumptions  were  made  in  SPAERS  regarding  the 
following  entities: 

Engines /Air era  ft 

SPAERS  assumes  that  engines  are  not  removed  from  one 
inoperative  aircraft  to  put  another  aircraft  up  (no  cannibali- 
zation). Only  serviceable  engines  removed  from  stock  may  be 
used  in  this  capacity.  Nor  is  any  distinction  made  between 
engines  in  stock  in  terms  of  their  life  or  age.  The  choice 
of  which  engine  services  which  aircraft  is  totally  random  in 
that  a failure  time  is  generated  only  after  installation. 
Engines  which  are  servicing  aircraft  with  same  combination 
codes  will  deteriorate  at  the  same  rate  (flying  hours  between 
removals).  Likewise,  two  engines  on  the  same  aircraft  will 
deteriorate  at  the  same  rate. 

The  engine  failure  (demand)  process  is  assumed  to  be 
a Poisson  process,  implying  exponential  failure  times.*  A 
Poisson  process  also  supports  other  basic  assumptions, 
namely,  that  engines  will  be  found  uniformly  distributed 
throughout  system  pipelines.  Secondly,  the  memoryless 
property  of  the  exponential  distribution  allows  failure  times 
to  be  generated  according  to  the  same  distribution  regardless 
of  the  engine's  past  experience. 


*NAVMET  argues  against  Poisson  demands  because  of 
the  uncertainty  in  predicting  the  demand  rate. 
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Repair  times  and  travel  times  between  bases  are 
assumed  to  be  constant.  This  serves  to  reduce  the  system 
variance  and  improve  the  statistical  accuracy  of  the 
performance  measures. 

Average  repair  times  were  interpreted  to  include  the 
waiting  time  which  an  engine  would  experience  in  the  repair 
function.  For  this  reason,  the  repair  function  is  assumed 
to  have  an  infinite  number  of  repair  facilities  and  pro- 
cessing times  will  include  average  waiting  times. 

The  discipline  for  filling  backorders  at  a base  is 
first  in  first  out.  At  bases  which  are  not  ports,  no  engine 
may  be  In  stock  while  a backorder  exists  at  that  base. 

Aircraft  flying  rates  are  given  in  flying  hours  per 
month  for  a given  aircraft  a given  base.  In  determining 
a failure  date,  SPAERS  assumes  this  flying  rate  to  be 
uniformly  distributed  throughout  the  month. 

Carriers 

The  repair  system  is  described  as  a tree  in  which 
all  failures  of  type  occurlng  at  base  i are  always 
repaired  at  the  same  base.  In  the  case  of  a port  which 
supports  no  flying  or  repair  activity,  the  same  relationship 
holds  as  if  the  carrier  failures  occurred  at  the  port.  That 
is  to  say,  that  .v  degree  failures  which  are  delivered  to 
the  port  from  any  carrier  are  sent  to  the  same  repair  base. 

The  fact  that  a carrier's  flying  activity  is  only 


active  during  deployment  supports  the  following  assumptions. 
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1)  No  failures  can  occur  during  a carrier's  port 

time . 

2)  Statistical  observations  are  made  on  a carrier 
only  during  deployment.  (Backorders  or  other  system  condi- 
tions do  not  affect  carrier's  performance  during  its  docking.) 

Although  a carrier's  flying  activity  is  dormant  in 
port,  its  repair  function  is  not.  It  is  not  unlikely,  then 
for  an  engine  to  become  available  on  board  a carrier  through 
the  efforts  of  the  on-board  repair  function.  All  other 
resupply,  however,  is  accomplished  by  the  port  concurrent 
with  the  time  of  deployment.  No  resupply,  external  to  that 
of  the  carrier's  own,  can  be  performed  during  its  deployment. 

Engines  on  board  a docked  carrier  cannot  be  used  by 
the  port  to  supply  an  out  going  carrier,  nor  can  it  use 
those  engines  for  its  own  flying  activity. 

SPAERS  Incorporates  some  variability  in  a carrier's 
port/deployment  schedule  by  using  the  normal  distribution. 

It  was  believed  that  a fixed,  non-varying  carrier  schedule 
would  make  the  demand  process  too  regular.  No  information 
is  readily  available  on  carrier  schedules  except,  perhaps, 
for  average  port  or  deployment  times.  So  the  intuitive  appeal 
in  using  the  normal  distribution  to  vary  the  schedule  lies 
in  the  ability  to  control  the  mean  and  variance  according  to 
the  information  available,  without  introducing  any  special 
assumptions  about  the  particulars  of  the  scheduling  activity.* 

*The  ability  to  use  fixed  schedules  is  available. 

(See  Appendix  B on  user  availability.) 
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SPAERS  event  scheduling  Introduces  some  assumptions 
concerning  the  repair  system  In  the  chance  that  two  events 
may  occur  simultaneously  In  the  simulation.  The  event 
priority  Is  as  follows: 

1)  Engine  becomes  available. 

2)  An  aircraft  falls. 

3)  A carrier  will  dock  or  deploy. 

4)  A sampling  event  occurs. 

Engine  availability  is  executed  first  to  avoid  the 
unnecessary  deactivation  of  an  aircraft  or  the  deployment  of 
an  under-stocked  carrier. 

An  aircraft  is  allowed  to  fail  before  a carrier 
docks  to  avoid  the  event  that  an  engine  fails  on  board  a 
docked  carrier. 

The  sampling  activity  is  considered  last  so  that  all 
pending  events  may  occur  before  the  system  statistics  are 
calculated . 

Simultaneity  of  events  is  a rare  event,  both  in 
SPAERS  and  in  reality,  but  this  priority  rationale  prevents 
any  difficulty  In  executing  the  simulation  events. 

B.  Initial  Conditions 

As  in  many  simulation  studies,  SPAERS  is  concerned 
with  Investigating  the  'long  run'  or  'equilibrium'  behavior 
of  the  system.  Unfortunately,  a simulation  model  must  be 
started  and  stopped.  Because  of  the  factitiousness  Introduced 
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by  the  abrupt  beginning  of  a simulation  run,  the  performance 
of  the  simulated  system  does  not  become  representative  of  the 
corresponding  real-world  system  until  it  has  essentially 
reached  a "steady-state"*  condition.  Thus,  we  must  be 
concerned  with  the  problem  of  how  to  obtain  data  that  are 
relevant  for  predicting  the  'steady-state'  behavior  of  the 
real  system.  Tn  other  words,  we  would  like  to  observe  data 
from  a system  operating  under  'chance  causes' — random 
variation,  and  not  a system  whose  underlying  probabilistic 
behavior  is  dependent  on  initial  conditions. 

There  are  two  traditional  ways  of  dealing  with  the 
'initial  conditions'  problem: 

• Run  the  simulation  model  for  some  time  without 
collecting  data  until  the  simulated  system  has 
"essentially"  reached  a state  representative  of 
steady-state  conditions.  This  time  is  usually 
referred  to  as  a 'run-in'  or  'stabilization'  period. 

• Start  system  in  a state  as  representative  of  steady- 
state  conditions  as  possible  to  minimize  the  required 
length  of  the  'run-in'  period. 

It  is  difficult  to  estimate  how  long  this  'run-in' 
period  needs  to  be.  Furthermore,  some  data  are  not  available 
to  fully  start  the  system  in  a state  representative  of 


*The  term  "steady-state"  has  the  probability  theory 
definition  of  describing  the  state  of  a system  as  having 
reached  a limiting  equilibrium  probabilistic  behavior, 
independent  of  time. 
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steady-state  conditions. 

The  authors  have  Included  in  SPAERS  a huerlstlc 
method  for  calculating  a 'run-in'  period.  Also,  some 
knowledge  is  known  about  the  system  to  accurately  initialize 
portions  of  the  system  state  sufficiently  close  to  steady- 
state  conditions.  They  are  discussed  below. 

Initial  States  of  Certain  Entities: 

• Aircraft  #1,  i = 1,  2,  • • *,  total  number  of  aircraft 
Since  no  data  are  available  on  how  many  aircraft  on 
the  average  are  down  at  any  particular  base,  all 
aircraft  are  initialized  as  being  "operative."  This 
means  the  down  aircraft  queue  is  empty  and  the 
falltime  queue  is  full.  Thus,  initially  SPAERS 
calculates  the  life  of  each  engine  on  aircraft  i.  The 
previous  history  of  an  engine  need  not  be  considered 
since  the  exponential  distribution  has  the  memoryless 
property . 

A failure  time  must  also  be  computed  for  aircraft  i. 

• Base  , 1=1,  2,  • • *,  total  number  of  bases. 

Since  all  aircraft  are  initialized  as  being  in  an 
'operative'  state,  the  number  of  backorders  and  the 


r number  of  down  aircraft  at  base  ^ are  both  equal  to 

zero.  Furthermore,  the  stock  level  at  base  ^ equals 
I the  number  of  spare  engines  belonging  to  base  1. 

There  are  two  components  of  base  1 stock: 

f 

: (1)  Safety  stock 

i. 


I 
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(2)  Pipeline  stock 

Both  of  these  quantities  are  either  computed  in  the 
initialization  using  the  *DOD  requirements'  method  or 
are  data  inputs  (user  option — see  Appendix  B).  The 
on-hand  stock  level  of  base  ^ is  initialized  to  equul 
the  safety  stock.  Each  engine  of  pipeline  stock,  on 
the  other  hand,  is  initialized  uniformly  in  base  ^'s 
pipeline.  Thus,  by  scheduling  the  availability  of 
these  engines  to  base  i in  the  future,  we  are  closer 
to  attainlrig  steady-state  conditions;  that  is,  in  the 
real  system  you  expect  to  observe  a certain  number  of 
engines  in  base  ^'s  pipeline  (resupply)  and  rarely 
v/ill  all  engines  exist  as  on-hand  stock. 

• Carrier  #1,  1 = 1,  2,  ' ' ’,  total  number  of  carriers 
The  state  of  carrier  ^ (deployment  or  docked)  is 
randomly  determined  with  Probability  {STATE(i)  = 

'AT  SEA’}  - mean  time  at  sea/(mean  time  at  sea  + mean 
time  at  port)  and  Probability  (STATE(i)  = 'AT  PORT'}  = 

1 - Probability  {STATE(l)  = 'AT  SEA'}. 

The  time  of  next  state  change  of  carrier  ^ is 
initialized  at  a fraction  of  its  first  schedule  time. 
This  schedule  time  may  be  deterministic  or  probabilistic 
(user  option — see  Appendix  B). 

Thus , 

Time  of  next  deployment  = U • first  deployment 

time 


if  STATE(i) 


'AT  PORT' 


l! 


[i 


i 
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or 

Time  of  next  docking  = U • first  docking  time, 
if  STATE(i)  = 'AT  SEA' 

where  U is  (0,  1)  uniformly  distributed  random 
variable . 

RUN-In  period  : 

At  any  point  in  time  an  engine  may  exist  as 
on-hand  stock,  pipeline  stock  or  as  an  operating  unit  on 
an  aircraft.  The  spare  engines  have  appropriately  been 
classified  and  located  in  the  system  accordingly  as  an 
attempt  to  achieve  steady-state  conditions.  Due  to  the 
lack  of  information,  we  were  forced  to  initialize  the 
remaining  engines  in  the  system  as  'operating  units  on 
aircrafts.'  To  alleviate  this  artificial  condition,  it 
seems  intuitively  appealing  to  allow  all  these  engines 
to  fail  at  least  once  before  any  collection  of  system 
data.  By  allowing  all  engines  in  the  system  to  exist  as 
on-hand  or  pipeline  stock  or  as  operating  units,  the 
system  has  come  one  step  closer  to  achieving  steady- 
state  conditions. 

Specifically,  the  run-in  period  is  computed  by: 

(1)  Calculating  the  average  time  it  takes  for  all 
'operating'  engines  to  fail  at  base 

(2)  Obtaining  the  maximum  time  over  all  'flying' 
bases . 


1 


(3)  Advancing  the  calculation  to  the  nearest  whole 
day . 


It  should  be  noted  that  during  SPAERS  initialization  a 
report  of  system  parameters  and  spare  stock  levels  is  printed. 

C . SPAERS  Options  and  Availabilities 

Variations  of  system  parameters  accomplished  through 
data  input  provide  the  user  the  opportunity  to  study  a 
system's  performance  and  sensitivity  under  different 
conditions . 

Some  examples  are: 

• Twisted  vs.  Continuous  Pipeline  Configurations 

• Probabilistic  vs.  Deterministic  Carrier  Deployment  and 
Dock  Times 

• If  probabilistic  carrier  times  are  employed,  then  the 
effect  of  high  variance  or  instability  may  be  studied. 

• The  effects  of  over  or  under  stocking  various  bases 
may  be  examined. 

• DODI-'4230 . ■4  [1]  calculated  spare  stock  levels  with  some 
confidence  level  vs.  user  data  input  spare  stock 
levels . 

• Effects  of  Including  (excluding)  a port  which  maintains 
neither  a flying  nor  a maintenance  activity. 

For  a detailed  discussion  oi.  these  and  other  SPAERS 


availabilities  the  reader  should  refer  to  Appendix  B. 
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D . Data  Output 


At  the  termination  of  the  simulation,  SPAERS  reports 
statistical  measures  which  summarize  base  performance 
activity  during  the  simulation  run. 

Measures  such  as  the  ready  rate,  fill  rate,  average 
on-hand  stock,  average  backorders,  average  down  aircraft, 
average  backorder  time,  average  down  aircraft  time  along 
with  non-zero  distribution  values  of  net  inventory  (on-hand 
stock — backorders)  and  down  aircraft  are  printed  for  each 
base.  Furthermore,  the  sample  mean  and  standard  deviation 
of  the  mean  along  with  a confidence  Interval  of  the  criterion 
(average  down  aircraft/total  number  of  aircraft)  for  each 
base  are  provided. 

Finally,  two  measures  of  system  criterion  are 
reported  to  the  user.  A detailed  discussion  of  the  meaning 
and  implications  of  these  measures  is  provided  in  Appendix  C. 


E.  SPAERS  System  Breakdown 
Entities  and  their  Attributes 


Entity : 
(i)  System — 


Attribute  : 

• Engine  Model  Type 

• Total  number  of  bases 

• Total  number  of  carriers 

• Total  number  of  aircraft 

• Total  number  of  aircraft  types 


j 

j 
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Entity : Attribute ; 

(11)  Engine  Model  Type — 

• Name 

• Total  number  of  'combination 

codes ' 

• Repair  time  for  each  degree  of 

failure 

(111)  Base  — 

• Name 

• Assigned  'base'  number 

• Assigned  'carrier'  number  ('O' 

If  ground  base) 

• On-hand  stock  level  at  any 

point  In  time 

• Number  of  backorders  at  any 

point  In  time 

• Total  number  of  backorders  up 

to  any  point  in  time 

• Number  of  down  aircraft  at  any 

point  In  time 

• Total  number  of  down  aircraft 

up  to  any  point  In  time 

• Number  of  aircraft  types 

• Number  of  aircraft  of  each  type 

• Total  number  of  aircraft 


Total  number  of  engines 


Entity : 


(iv)  Carrier — 


(v)  Aircraft — 


Attribute : 

• Travel  time  to  every  base  in 

system 

• Base  at  which  an  engine  will 

be  repaired  when  failure 
degree  i occurs.  (-1  ^ i - 

• Desired  criterion 

• Weight  o’f  importance 

• Highest  repair  degree 

• Time  of  next  deployment  or 

docking 

• State  (deployment  or  docking) 

at  any  point  in  time 

• Port  at  which  carrier  docks 

• Future  Deployment  and  Docking 

Schedules 

• Number  of  engines  that  have 

been  sent  off  to  be  repaired 
but  have  not  been  satisfied 

• Number  of  demands  (engines) 

that  were  not  satisfied  by 
carrier's  port  during  last 
'docking ' 


Assigned  number 


Entity ; 


Attribute : 


• Base  to  wnlch  aircraft  belongs 

• Type  of  aircraft 

• Next  failure  time 

• Remaining  life  time  of  engine ( s ) 

on  aircraft  at  time  of 
failure  (In  terms  of  flying 
hours ) 

• Flying  Activity  In  terms  of 

flying  hours  per  month 

(vl)  Type  of  Aircraft  — 

• Name 

• Assigned  number 

• Combination  code 

• Engine  capacity 

(vll)  Combination  Code — 

• Actual  code  number  used  by  NAVY 

• Assigned  number 

• Probability  of  degree 

failure  for  all  ^ (-1  ^1-3) 

• Mean  time  between  removal  (In 

terms  of  flying  hours) 

Events  and  their  Scheduling 

An  endogenous  event  takes  place  as  a result  of  Its 
being  scheduled  during  the  course  of  the  simulation 
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run.*  SPAERS  endogenous  events  are  as  follows: 

• The  time  at  which  an  engine  failure  causes  an  aircraft 
to  become  Inopei’atlve . This  event  is  appropriately 
referred  to  as  PLANEFAILEVENT . 

• The  time  at  which  an  engine  in  a particular  pipeline 
becomes  available  as  on-hand  serviceable  stock  at  a 
base.  This  event  is  referred  to  as  ENGAVAILEVENT . 

• The  time  at  which  the  state  of  a particular  carrier 
changes.  This  event  is  referred  to  as  SHIPEVENT. 

An  exogenous  event  is  made  to  occur  by  specifying  the 
type  of  event  and  simulated  time  at  which  it  is  to  occur 
through  data  input  prior  to  the  start  of  a simulation  run.* 
This  type  of  event  usually  serves  as  a means  of  controlling 
a deslx’ed  aspect  of  the  simulation  run.  SPAERS  exogenous 
events  are  as  follows : 

• The  time  at  which  appropriate  statistics  are  to  be 
gathered  for  analysis  purposes  (see  Appendix  AIII  for 
discussion).  This  event  is  referred  to  as  SAMPLE. 

• The  time  at  which  the  simulation  run  is  to  be 
terminated  and  final  statistics  are  reported.  This 
event  is  appropriately  labeled  REPORT. 

SPAERS  synchronization,  that  is,  the  identification 
of  the  event  occurring  next  and  its  time  of  occurrence,  is 
based  on  asynchronous  timing — the  search  for  a minimum  among 
event  times  (see  Flowchart  in  Appendix  All). 


*Reference  [6],  pp . 232-233. 
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To  illustrate  SPAERS'  treatment  of  the  two  pipeline 
configurations,  an  actual  repair  system  was  analyzed.  The 
U.  S.  Navy's  PASER  data  were  used  to  calculate  DODI  spare 
stock  levels.  PASER  also  provided  system  parameters  that 
served  as  SPAERS  input. 

AIRPAC  Command's  Engine  type  model  T58GE-series 
8B  and  8F  in  the  North  Island  Naval  Aircraft  Rework  Facility 
System  (Combination  Codes  731j  733  and  735)  defined  the 
repair  system  which  was  simulated  by  SPAERS  in  this  analysis. 

The  DODI  calculation  was  performed  using  the 
following  consideration.  According  to  PASER,  ports  which  do 
not  fly  or  repair  will  not  enter  into  spare  stock  level 
computations.  That  is,  carriers  are  connected  directly  to 
their  repair  echelons,  ignoring  the  port.  The  implication 
here  is  that  the  port  is  not  equipped  to  satisfy  carrier 
demands,  but  the  carrier  repair  echelons  are. 

A 0.9  ready  rate  criterion  was  used  for  all  bases. 
SPAERS  showed  that  the  configuration  assumed  by  DODI  and 
stocked  accordingly,  did,  in  fact  produce  ready  rates  close 
to  0.9. 
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To  demonstrate  the  effect  of  Introducing  the  port  to 
a continuous  pipeline  situation,  the  following  configurations 
were  run  by  SPAERS. 

1)  Continuous  carrier  pipelines  without  ports 

(DOD  spare.s  calculation  included) 

2)  Continuous  carrier  pipelines  with  ports 

3)  Twisted  carrier  pipelines  with  ports 

Comparisons  of  the  down  aircraft  criterion  were  made 

in  each  case  using  the  test  of  the  hypothesis  that  the  means 
of  two  Normal  Distributions  are  equal,  assuming  that  the 
standard  deviations  are  not  known  and  not  necessarily  equal* 
(.05  level).  The  test  was  performed  between  corresponding 
bases  in  two  different  configurations. 

The  first  comparison  Involved  1)  and  2)  to  show  any 
significant  effect  on  performance  by  the  introduction  of 
ports  to  the  original  DODI  system. 

As  the  results  in  table  1 indicate,  only  Cubi  Point, 
a port,  shows  a significant  difference  in  performance.  In 
reality,  Cubi  Point  repairs  and  supports  its  own  flying 
activity.  By  virtue  of  its  port  activity,  it  experiences 
higher  engine  availability  by  storing  engines  for  carriers. 
But  these  engines  may  be  used  for  Cubi  Point's  own  flying 
activity  and  serve  to  decrease  its  down  aircraft  time. 

Another  comparison  was  made  between  configuration 


*Reference  [5],  pp.  240-2^12. 
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2)  and  3).  Here  the  singular  effect  of  a twisted  pipeline 
is  contrasted  with  the  continuous  pipeline.  Again,  Cubi 
Point  is  significantly  different  from  the  continuous 
situation.  The  effect  of  increased  engine  availability  at 
this  port  is  enhanced  by  the  twisted  pipeline  behavior  of 
carriers.  Since  carriers  demand  engines  before  they  deploy 
Instead  of  continuously  as  in  the  case  of  the  continuous 
pipeline,  the  engines  which  are  stocked  at  the  port  in 
anticipation  of  carrier  demands  may  also  be  used  for  the 
port's  flying  activity.  The  results  of  this  comparison  appear 
in  table  2. 

Table  2 also  shows  that  the  CV-PAC-1  carrier 
performs  significantly  worse  in  the  twisted  pipeline.  It 
was  noted  that  this  carrier  supports  a more  extensive  flying 
activity  with  a yearly  failure  rate  of  28.17  engines 
than  other  system  carriers.  These  other  carriers  produce 
3 or  4 failures  per  year. 

The  statistical  comparison  between  1)  and  3) 
compounds  the  situation  of  ports  and  twisted  carrier 
pipelines.  Once  again,  the  port  and  its  most  active  carrier 
demonstrate  significantly  different  performances  (table  3). 

Another  system.  Identified  by  engine  model  type 
J52P  - series  8A  and  8B  in  the  Jacksonville  Naval  Aircraft 
Facility  System  (Airlaut  Command,  Combination  Codes  102  and 
104)  was  observed  in  the  same  manner.  Stock  levels  close  to 
DOD  requirements  were  used  for  each  base  and  in  both  the 
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continuous  and  the  twisted  pipeline  situations  were  run 
(including  ports).  The  three  carriers  in  this  system  were 
particularly  active,  with  51.6  yearly  engine  failures.  As 
table  4 shows,  carriers  are  significantly  different  in  terms 
of  their  performance  measures  in  each  configuration. 

It  should  be  noted  that  no  other  repair  system 
locations  exhibit  criterion  differences  using  this  test  of 
hypothesis  for  equal  means  in  either  of  the  two  systems 
observed . 
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TABLE  1 

Comparison  of  Mean  Criterion — Continuous  Pipeline 
Without  Ports  vs.  Continuous  Pipeline  with  Ports: 
Engine  Model  Type  T56GE  83,  8F 


mean  criterion  in  Continuous 
Pipeline  Configuration  without 
non-f lying/maintenance  ports 

standard  deviation  of 

mean  criterion  in  Continuous  Pipe- 
line Configuration  with  non-flying/ 
maintenance  ports 

standard  deviation  of  x- 


Spare  Stock 
Level 


Test  Signifi- 
Statlstlc  cant? 


Adak 

1 

0.04428 

0.0051 

0.04197 

0.00581 

0.298 

Alameda 

0.01397 

0.00216 

0.01179 

0.00216 

0.713 

Barbers 

Point 

2 

0.01145 

0.00243 

0.00685 

0.00134 

1.658 

Cubl 

Point 

4 

0.01314 

0.00414 

0.02821 

0.0053 

2.24 

CV-PAC-1“ 

3 

0.01211 

0.00236 

0.01133 

0.00331 

0.1918 

CV-PAC-il* 

1 

0.006 

0.00353 

0.00315 

0.00171 

0.7266 

Imperial 

Beach 

0 

0.01163 

0.00494 

0.01559 

0.0036 

0.6478 

Kadena- 

AFB 

0 

0.03174 

0.00924 

0.0233 

0.00703 

0.7269 

Kaual- 

PMRF«* 

NAV-PAC-1* 

1 

0.01931 

0.00968 

0.00587 

0.0044 

1.264 

NAV-PAC-2* 

1 

0.02074 

0.01112 

0.00652 

0.00533 

1.153 

NAV-PAC-3* 

1 

O.OI6O9 

0.00791 

0.01117 

0.00554 

0.509 

NAV-PAC-4* 

1 

0.00815 

0.00566 

0.00650 

0.00453 

0.227 

NAV-PAC-5* 

1 

0.01894 

0.009 

0.01156 

0.00546 

0.701 

North 

Island 

10 

0.00473 

0.00071 

0.00494 

0.00078 

0.199 

Point  Hugu 

1 

0.00328 

0.00177 

0.00289 

0.00149 

0.173 

Wldfey 

Island 

1 

0.00771 

0.00248 

0.00436 

0.00149 

1.157 

Iwakunl-MC 

1 

0.01781 

0.0051 

0.02116 

0.00447 

0.4939 

Kaneohe-MC 

1 

0.00248 

0.00105 

0.00327 

0.00115 

0.507 

* Carrier 

Aircraft  do  not  ’fly* 


TABLE  2 


Comparison  of  Mean  Criterion — Continuous  Pipeline 
With  Ports  vs.  Twisted  Pipeline  With  Ports: 
Engine  Model  Type  T58GE  8b,  8F 


Z,  : mean  criterion  in  Continuous  Pipeline 

Configuration  with  non-f lying/mainten- 
ance ports 

oT^ : standard  deviation  of 

X"p : mean  criterion  in  Twisted  Pipeline 

Configuration  with  non-flying/ 
maintenance  ports 

Op'-  standard  deviation  of  Xp 


Base  Spare  Stock 

Level 

: 

X2 

^1 

Op  Test 

Statistic 

Signif  i' 
cant? 

Adak 

1 

0.04197 

0.04707 

0.00581 

0.00672 

0.5741 

Alameda 

0.01129 

0.01008 

0.00216 

0.00208 

0.4035 

Barbers 

Point 

2 

0.00685 

0.00770 

0.00134 

0.00165 

0.39989 

Cubl 

Point 

0.02621 

0.00169 

0.00530 

0.00073 

4.9569 

/ 

CV-PAC-1* 

3 

0.01133 

0.03818 

0.00331 

0.00551 

4.177 

/ 

CV-PAC-i|* 

1 

0.00315 

0 . 00847 

0.00171 

0.0034 

1.3978 

Imperial 

Beach 

0 

0.01559 

0.00890 

0.0036 

0.00353 

1.326 

Kadena-AFB 

0 

0.02330 

0.02111 

0.00703 

0.00655 

0.2279 

Kaual- 

PMRF*» 

does  not 

fly 

NAV-PAC-1* 

1 

0.00587 

0.02187 

0.0044 

0.00856 

1.662 

NAV-PAC-2* 

1 

0.00652 

0.01802 

0.00533 

0.00594 

1.441 

NAV-PAC-3* 

1 

0.01117 

0.00460 

0.00554 

0.00342 

1.009 

NAV-PAC-^* 

1 

0.0065 

0.01122 

0.00453 

0.00551 

0.661 

NAV-PAC-5* 

1 

0.01156 

0.02503 

0.00546 

0.01214 

1.012 

North 

Island 

10 

0.00494 

0.00601 

0.00078 

0.00069 

1.027 

Point  Mugu 

1 

0.00289 

0.00310 

0.00149 

0.00146 

0.10067 

Wldfey 

Island 

1 

0.00436 

0.00435 

0.00149 

0.00118 

0.00526 

Iwakuni-MC 

1 

0.02116 

0.02605 

0.00447 

0.00575 

0.6714 

Kaneohe-MC 

1 

0.00327 

0.00703 

0.00115 

0.00209 

1.576 

* Carrier 

**  Aircraft  do  not  'fly' 


TABLE  3 


Comparison  of  Mean  Criterion — Continuous  Pipeline 
Without  Ports  vs.  Twisted  Pipeline  With  Ports: 
Engine  Model  Type  T58GE  8B,  8F 


X, : mean  criterion  in  Continuous  Pipeline 

Configuration  without  non-flying/ 
maintenance  ports 

: standard  deviation  of 

X2’.  mean  criterion  in  Twisted  Pipeline 
Configuration  with  non-flying/ 
maintenance  ports 

02".  standard  deviation  of  X'2 


Base  Spare  Stock  x. 
Level 


Xp  a.  Op  Test 

Statistic 


Adak 

1 

0.0i<il28 

0.04707 

0.0051 

0.00672 

0.33 

Alameda 

4 

0.01397 

0.01008 

0.00216 

0.00208 

1.297 

Barbers 

Point 

2 

0.01145 

0.00770 

0.00243 

0.00165 

1.276 

Cubi 

Point 

H 

0.01314 

0.00169 

0.00414 

0.00073 

2.72 

CV-PAC-1* ** 

3 

0.01211 

0.03818 

0.00236 

0.00551 

4.349 

CV-PAC-4* 

1 

0.006 

0.00847 

0.00353 

0.0034 

0.503 

Imperial 

Beach 

0 

0.01163 

0.00890 

0.00494 

0.00353 

0.449 

Kadena- 

AFB 

0 

0.03174 

0.02111 

0.00924 

0.00655 

0.919 

Kaual- 

PMRF*» 

does  not 

fly 

NAV-PAC-1* 

1 

0.01931 

0.02187 

0.00968 

0.00856 

0.198 

NAV-PAC-2* 

1 

0.02074 

0.01802 

0.01112 

0.00594 

0.216 

NAV-PAC-3* 

1 

0.01609 

0.0046 

0.00791 

0.00342 

1.333 

NAV- PAC-H* 

1 

0.00815 

0.01122 

0.00566 

0.00551 

0.388 

NAV-PAC-5* 

1 

0.0189'* 

0.02503 

0.009 

0.01214 

0.4029 

North 

Island 

10 

0.00473 

0.00601 

0.00071 

0.00069 

1.293 

Point  Mugu 

1 

0.00328 

0.0031 

0.00177 

0.00146 

0.078 

Wldf ey 
Island 

1 

0.00771 

0.00435 

0.00248 

0.00118 

1.22 

Iwakuni-MC 

1 

0.01781 

0.02605 

0.00510 

0.00575 

1.072 

Kaneohe-MC 

1 

0.00248 

0.00703 

0.00105 

0.00209 

1.945 

* Carrier 

**  Aircraft  do  not  'fly' 
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TABLE  4 


Comparison  of  Mean  Criterion — Continuous  Pipeline 


With  Ports  vs 
Engine 

. Twisted  Pipeline  With  Ports: 
Model  Type  J52P  8a,  8B 

Xp: 

mean  criterion  in  Continuous  Pipeline 
Configuration  with  non-flying/ 
maintenance  ports 

Op: 

standard  deviation  of  x^ 

: 

mean  criterion  in  Twisted  Pipeline 
Configuration  with  non-flying/ 
maintenance  ports 

02 : 

standard  deviation  of  X2 

1 

I 

I 

i 


] 

I 

J 

D 
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Base  Spare  Stock 

Level 

; x^ 

X2 

^1 

^2  , 

Test 

Statistic 

Slgnifl 

cant? 

CV-LANT  1* 

5 

0.00357 

0.01909 

0.00109 

0.00206 

-6.6592 

/ 

CV-LANT  2* 

5 

0.00^77 

n. 01535 

0.00081 

0.00159 

-5.93 

/ 

CV-LANT  3* 

0.00750 

0.01947 

0.00124 

0.00204 

-5.014 

/ 

Cherry 

Point 

6 

0.00897 

0.00838 

0.001 

0.00078 

0.4652 

South 

Weymouth 

1 

0.01809 

0.01840 

0.00153 

0.00171 

-0.135 

Oceana 

25 

0.02737 

0.02624 

0.00214 

0.00178 

0.4059 

Key  West 

1 

0.00915 

0.00893 

0.00272 

0.00267 

0.0577 

Norfolk 

2 

0.01326 

0.01360 

0.00226 

0.00232 

-0.1049 

Paxtuxent 

Riv 

1 

0.00116 

0.00323 

0.0071 

0.00183 

-0.1895 

Willov; 

Grove 

1 

0.02324 

0.02118 

0.00231 

0.00204 

0.668 

Memphis 

1 

0.0036 

0.00644 

0.00157 

0.00241 

-0.428 

* Carrier 
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RECOMMENDATIONS  AND  GENERAL  REMARKS 

The  previous  section  demonstrated  the  use  of  SPAERS  in 
computing  performance  of  repair  system  configurations. 

Based  on  the  runs  of  two  actual  repair  system,  it  was 
observed  that  some  carriers  and  ports  in  a twisted  pipeline 
perform  significantly  differently  from  the  sane  carriers  and 
ports  in  the  continuous  case.  The  remainder  of  the  system  (in 
these  cases)  appeared  to  perform  independently  of  carrier 
behavior. 

It  was  also  suggested  that  the  extent  to  which  a carrier 
in  the  twisted  pipeline  configuration  will  perform  differently 
from  the  continuous  one  is  somehow  indicated  by  the  carrier's 
failure  rate.  The  nature  of  the  twisted  pipeline  carrier 
activity  resembles  that  of  a periodic  review  Inventory  system 
where  each  inventory  review  occurs  during  a carrier's  dock  time. 
When  this  is  contrasted  with  the  continuous  review  analog  for 
the  continuous  pipeline,  inventory  requirements  for  the  twisted 
pipeline  would  be  higher.  If  the  carrier's  flying  activity  (and 
so  the  failure  rate)  is  low,  the  two  inventory  systems  have 
similar  stock  requirements . On  the  other  hand,  if  activity  is 
high,  the  periodic  review  will  require  more  stock  to  adequately 
support  its  flying. 
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Stronger  statements  about  expected  performances  for  a 
•'i  general  system  would  require  more  extensive  investigation. 

The  authors  recommend  the  following: 
i*  1)  A variety  of  systems  should  be  compared  in  this  way  to 

support  some  general  statements  concerning  carrier  activity. 

DOD  stock  levels  are  adequate  for  this  investigation, 
j 2)  NAVMET  levels  should  be  used  as  input  stock  levels  to 

i-  i 

SPAERS  simulation  of  repair  systems,  noting  the  difference 
bebween  this  system  performance  with  DOD  levels. 

3)  Ports  and  their  activity  with  carrier  operations  should 

I 

be  considered  in  analytical  spare  stock  level  allocation  to 
account  for  their  role  in  carrier  performance. 

I 1 
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Appendix  A:  SPAERS  Delineation 
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AI.  Linked  Listing/Queue  Mechanism 
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Linked  listing  is  a common  technique  used  in  SPAERS 
even  scheduling.  The  invariant  relation  in  a linked  list  is 
the  fact  that  all  the  entries  in  the  list  are  ordered 
according  to  some  discipline.  The  method  of  ordering  relies 
on  a system  of  "pointers,"  or  attributes  of  each  list  entry 
which  indicates  the  address  of  the  next  entry  in  the  list. 

Each  entry  is  identified  by  only  one  address.  Its  position 
in  the  list  is  identified  by  the  "pointers"  or  system  links. 

The  application  of  linked  lists  to  computer 
Information  management  is  crucial  to  efficiency.  A new  entry, 
for  instance,  does  not  require  re-shlftlng  any  part  of  the 
list.  Only  two  "pointers"  need  to  be  changed.  Similarly, 
removal  of  an  item  from  the  list  requires  only  "pointer" 
changes.  The  expected  number  of  comparisons  in  a list  search 
is  the  same  as  any  other  linear  search  (1/2  N,  where  N is  the 
average  size  of  the  list). 

Let  the  following  illustration  represent  computer 
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storage  In  the  form  of  an  array. 

For  the  purpose  of  this  example,  ordering  attributes 
can  be  event  times  which  will  enter  and  leave  the  list. 

To  maintain  the  list  ordering,  an  auxiliary  array  is 
required.  Each  entry  in  this  array  indicates  the  address  of 
the  entry  following  it  in  the  list: 


Scheduler 


0 

1 

2 

3 

4 

99  9 9 

11:00 

11:30 



11:20 

11:25 

< — address 


* att  ri but  e 

(time) 


0 

1 

r 

2 

3 

4 

1 

3 

0 

4 

2 

EMPTY  ADDRESS--5 

act i ve 


-address 


c — attri  bute 
(address  of 
attribute 
immediate! 


the 

y 


empty 


follow  i ng ) 


Note  that  the  linked  list  may  be  divided  into  an 
active  and  an  empty  section.  The  active  section  consists  of 
those  addresses  whose  attributes  are  in  the  schedule.  The 
empty  section  consists  of  those  addresses  whose  attributes 
are  not  part  of  the  active  list.  Both  sections  are  linked 
and  the  leading  element  of  both  sections  can  be  identified. 

It  is  convenient  to  have  the  zero  entry  in  the  link 
system  so  that  its  link  attribute  Identifies  the  address  of 
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of  the  leading  entry  in  the  active  list.  The  link  of  the 
final  entry  in  the  schedule  is  back  to  the  zero  address,  thus 
forming  a linked  cycle.  The  zero  address  schedule  attribute 
should  be  of  a magnitude  so  that  it  will  always  be  at  the 
end  of  the  schedule.  In  this  way,  the  zero  entry  acts  as  a 
schedule  "anchor"  in  that  any  other  entry  in  the  list  must 
proceed  it,  and  its  link  attribute  always  identifies  the  top 
of  the  list.  (A  special  case  worth  mentioning  occurs  when 
the  zero  link  attribute  identifies  its  own  address,  signify- 
ing an  empty  schedule  list.) 

To  add  an  element  to  a linked  list : 

1)  enter  the  attribute  in  the  first  address  available 
in  the  empty  list,  retaining  the  address. 

2)  change  the  empty  address  identifier  to  the  next  empty 
address . 

3)  Identify  the  point  of  insertion  for  the  new  item  by 
searching  the  ordering  attribute  list  and  utilizing 
the  link  system.  The  address  of  the  proceeding 
entry  will  identify  an  insertion  point  (B  in  the 
figure ) . 

4)  link  the  new  item  with  the  old  address  link  of  B. 

Link  B to  the  address  of  the  new  item.  Insertion  is 
complete  and  the  invariant  relation  still  holds. 


AM 


—old  link 

new  link 


etween  B and  A 

In  deleting  an  element  N from  the  list, 

1)  Locate  the  Item  to  be  deleted  by  identifying  its 
address  N and  the  address  of  Its  preceedlng  link  B. 
(In  event  schedules,  the  outgoing  Item  is  usually 
the  leading  entry.) 

2)  Link  B to  A retaining  the  address  of  N. 

3)  Link  N with  the  address  of  the  first  empty  element 
and  Identify  N as  the  new  first  empty  element. 

Variations  on  this  technique  Involve  the  use  of 
"backpointers"  to  facilitate  a search  which  starts  at  the 
end  of  the  list.  Other  attribute  lists  may  be  used  to  de- 
scribe the  elements  which  are  found  in  a specific  ordering. 

SPAERS  event  scheduling  relies  on  linked  list  systems 
(later  referred  to  as  queues).  Engine  availability,  aircraft 
failure,  and  carrier  events  are  each  queued  in  their  own  link 
systems  in  order  of  their  schedule  occurrence.  The  leading 
entry  in  each  link  system  is  the  next  event  which  will  occur 
of  that  type.  In  order  to  determine  the  next  event,  only  the 
leading  elements  of  each  link  system  require  comparison. 


To  insert  N 
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Aircraft  failure  queue 

Any  aircraft  which  is  not  in  the  down  state  will  be 
assigned  a failure  date.  All  operational  aircraft  can  be 
found  in  the  failure  queue  and  ordered  by  failure  dates. 

Most  of  the  system  aircraft  are  expected  to  be  in  an 
operational  state  at  any  point  in  time  making  the  active 
portion  of  this  queue  large.  It  can  be  seen  that  by 
permanently  addressing  the  aircraft's  failure  time  by  its 
actual  aircraft  number,  maintenance  of  an  empty  queue  can  be 
avoided . 


-a/c  number 
-tail  date 


4 a/c  with  next 

fail  date 


Note  that  in  the  above  queue  link  system,  aircraft 
numbers  "2"  and  "6"  are  in  the  down  state  and  are  not  part 
of  the  link  system. 

Identifiers  used  in  the  failure  queue. 

ACFAILQFPT  (N)  is  the  address  (a/c  number) 
which  will  fail  after  aircraft  N. 

ACFAILQBPT  (N)  is  the  address  of  the  aircraft 
which  will  fail  before  N. 

ACFAILQBOT  identifies  the  aircraft  number  to 
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fall  last. 

BASEQ  Identifies  the  base  location  of  an  a/c . 


Availability  Queue 


Availability  Queue  sequences  the  times  of  engine 
availability  throughout  the  system.  Since  engines  are  not 
specifically  identified  by  number,  they  have  no  readily  made 
address  as  in  the  case  of  the  failure  queue.  An  engine  which 
enters  the  system  pipeline  as  a result  of  a demand  or  a 
failure  is  queued  in  order  of  its  availability  date.  The 
base  at  which  an  engine  will  become  available  is  also  an 
attribute  associated  with  the  address  of  a queued  engine. 

Relevant  identifiers  for  this  system  include: 

BASENGAVAIL  (N) — the  base  where  the  engine 
identified  by  address  N will  become  available. 

ENGAVAILBPT  (N) — the  address  of  the  engine  which 
will  become  available  in  the  system  prior  to  the  engine  with 
address  N. 

ENGAVAILQFPT — the  address  of  the  engine  which 
will  become  available  in  the  system  after  the  engine  whose 
address  is  N. 

ENGAVAILQNT — identifies  the  address  of  the  first 
available  array  position. 

ENGAVAILQHOLD — retains  addresses  during 
insertion  cr  extraction. 


, . 


ENGAVAILQBOT — identifies  the  engine  address 
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which  will  experience  the  last  availability  date. 

ENGAVAILTIME  (N) — availability  date  of  the 
engine  with  address  N. 

Ship  queue 

Ship  queue  orders  carrier  events  by  date.  All 
carriers  will  always  occupy  this  queue,  since  they  always 
have  an  event  trending.  In  this  case,  the  carrier  number  is 
its  queue  address. 

Identifiers  used  in  queue  maintenance  are: 

SHIPTIME  (N)  is  the  time  when  the  next  event 
(dock  or  deploy)  will  occur  for  carrier  N. 

SHIPQPT  (N)  Identifies  the  carrier  which  will 
experience  an  event  after  carrier  N. 

INSRT  and  LINSRT  retain  addresses  during 
reslnsertlon  process. 

Note  that  SHIPQPT  (0)  identifies  the  carrier 
with  the  earliest  future  carrier  event. 


I 

1 

1 
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Backorders 

It  was  stated  earlier  that  demands  that  cannot  be 
levelled  (filled)  immediately  are  logged  as  backorders  to  be 

filled  on  a first  come  first  served  basis  as  engine  avail- 

i 

[ ability  permits.  This  requires  that  the  backorder  commit- 

I ments  be  queued  at  each  base  in  the  order  of  v;hich  they  occur. 
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The  queued  attribute  is  that  base  to  which  the  backorder  is 
due . 

In  this  case,  queue  insertion  always  occurs  at  the 
"bottom"  of  the  list  and  removal  is  done  at  the  "top."  Thus 
only  the  leading  and  the  final  backordered  base  address  need 
to  be  identified  with  the  base  which  logged  the  backorder. 
The  mechanism  may  be  pictured  as  follows: 


BOQTOP  (N)  Identifies  joints  to)  the  address  of 
Base  N's  first  backorder  commitment. 

BOQBOT  (N)  identifies  the  address  of  the  base  N's 
last  backorder  commitment.  (Note:  If  both  the  above 

variables  are  zero  for  some  N,  no  backorders  are  outstanding 
at  that  base.) 

A special  case  of  the  backorder  situation  occurs 
when  a base  must  backorder  to  itself  (its  own  flying  activity). 

A self-backorder  is  the  same  as  a downed  aircraft,  and  is 
queued  as  both.  Similar  to  the  backorder  case,  down  aircraft 


r 
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^ are  re-installed  In  the  same  order  in  which  they  occurred. 

The  queuing  mechanism  has  a similar  operation  to  the  back- 

. order  mechanlfm  (see  page  Ad).  The  difference  with  the  down 

aircraft  queue  is  that  aircraft  addresses  are  actual  aircraft 
numbers  (as  in  the  case  of  the  aircraft  failure  event 

I schedule . 

[ DACQPT  (N)  indicates  the  next  aircraft  in  some 

5 

f base's  down  queue. 

i 

DACQTOP  (N)  identifies  the  first  down  aircraft 

[ number  at  base  N (zero  if  no  down  aircraft  exist  at  base  N). 

[ DACQBOT  (N)  Identifies  the  last  down  aircraft 

number  at  base  N (zero  if  no  down  aircraft  exist  at  base  N ) . 

i .'i 


All.  SPAERS  Flowchart 


The  task  of  developing  a computer  program  is 
facilitated  by  preparing  a layout  of  the  logical  flow  of 
operations.  This  layout  serves  as  a basis  for  which  the 
source  instructions  are  written. 

An  assortment  of  symbols  for  the  basic  types  of 
internal  operations  are  used  for  the  preparation  of  flowcharts. 

SPAERS  flowchart  symbols  are  as  follows: 

1.  The  slot: 


■> 


Used  to  denote  the  beginning  or  end  of  the 
program.  It  serves  as  a point  of  reference  and  does 
not  represent  any  particular  operation. 

2.  The  rectangle: 


Y=  4/X 


— > 


AlO 


All 


Used  to  Indicate  an  Internal  operation  which 
corresponds  to  one  or  more  assignment  statements. 
The  contents  of  the  box  may  be  either  of  symbolic/ 
mathematical  form  or  of  narrative  form. 

3.  A hexagonal  shaped  box: 


Used  to  describe  an  Internal  operation  requiring 
a decision.  The  decision  point  Is  stated  as  a 
yes-or-no  question. 


Ui 


BO  queue 


SHIP;  carrier  tfiat  is  changing 
states. 

NS(i):  nurrher  of  er>gint-s  that 
h .ve  been  sen t off  for 
repair  frorn  carrier  'i* 


fdIHSISN 


4 


1 

L 

T~ 

% * 

AIII.  Description  of  Procedures 
(1)  PLANEFAILEVENT 

This  procedure  is  Invoked  at  the  time  an  engine  falls 
which  renders  an  aircraft  inoperative. 

• ‘ Nomenclature ; 

• BASEFAILEDAC — the  base  number  of  the  base  at  which 
aircraft  has  failed. 

• DESTINBASE — the  base  number  of  the  base  at  which 

1 

failed  engine  is  to  be  repaired. 

• FAILEDAC — aircraft  number  of  failed  aircraft. 

• FAIL.DEGREE — failure  degree  of  engine 

• HD(I)--the  highest  degree  failure  which  base  I can 
repair.  (’4'  denotes  no  maintenance  activity  at 
base  ) 

• LIFENG(I,J) — the  number  of  flying  hours  before 

t h 

failure  of  an  engine  in  the  position  on  aircraft 
J.  (I  = 1,  2,  • • • , N:  N H engine  capacity.) 

• NEXTECH( I , J) --base  number  of  base  that  repairs  base 
^'s  degree  J failed  engines. 

• REPAIRTIME ( I ) — repair  time,  in  days,  of  an  I^^^  degree 
failed  engine. 

• SHIPNUM(I)--the  carrier  (ship)  number  of  base  I^  ('O' 
if  a ground  base). 

Al6 
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• STOCKLEVEL( I ) — number  of  serviceable  engines  on  hand  • ; 

i 

at  base  I.  i 

- \ 

• CCDISTdjJ) — cumulative  distribution  of  failure 
degree  of  combination  code  I_,  i.e.,  Pr  {FAILDEGREE 
^ J}  for  J = -1,  0,  1,  2,  3. 

• DUMMYTIME — time  remaining  in  carrier’s  deployment  if 
BASEFAILEDAC  is  a carrier  in  'Twisted  Pipeline 
Configuration. ' 

— 0 otherwise. 

Synopsis ; 

—Determine  FAILEDAC  and  BASEFAILEDAC  from  FAILTIME 
and  BASEQ  queues,  respectively. 

— Remove  FAILEDAC  from  PAILTIME  queue  by  appropriately 
rearranging  the  pointer  mechanism. 

— Determine  FAILDEGREE: 

• Generate  a random  variable,  U,  uniformly 
distributed  between  (0,  1). 

• FAILDEGREE  will  then  equal  the  largest  Integer 
K,  such  that  {K;  U ^ CCDIST(I,K), 

K = -1,  0,  1,  2,  3>.  IE  combination  code  of 
the  aircraft  type  of  FAILEDAC. 

— Attempt  to  make  FAILEDAC  operational. 

• If  STOCKLEVEL  (BASEFAILEDAC)  > 0,  then; 

1)  Remove  a repaired  engine  from  stock  and 

place  the  engine  on  FAILEDAC. 


2)  Determine  the  life  of  the  engine,  LIFENG 
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*•  (1,  FAILEDAC)  (l.e.,  number  of  flying 

hours  before  failure),  by  Invoking  the 

* ^ 

RNG  function  procedure  (see  (x)  below). 

^ 3)  Determine  the  next  failure  time  (In  days) 

of  the  aircraft  by  Invoking  PAILGEN 
procedure  with  parameters:  (FAILEDAC, 

SHIPNUM(BASEFAILEDAC) , BASEFAILEDAC ) 

(see  (vll)  below). 

• If  STOCKLEVEL( BASEFAILEDAC)  = 0,  then: 

1)  Place  FAILEDAC  In  the  down  aircraft 
queue  utilizing  DACQPT,  DACBOT  and  DACTOP. 

2)  Record  a backorder  to  BASEFAILEDAC  from 
BASEFAILEDAC  by  utilizing  BASEBO  queue. 

3)  Queue  backorder  accordingly  by  Invoking 
BACKORDQ  procedure  with  parameter: 

(BASEFAILEDAC)  (see  (lx)  below). 

— Engine  Transaction 

• If  HD(DESTINBASE)  > FAILDEGREE,  l.e., 

DESTINBASE  cannot  repair  the  engine  as  in  the 
case  of  a port  which  has  neither  maintenance 
nor  flying  activity  but  Is  required  to  serve 
as  an  intermediary  point  for  all  Its  carriers 
(user  option — see  Appendix  B) , then  provisions 
must  be  made  for  demand  accounting: 


1)  The  transaction  of  the  failed  engine 
will  be  between  the  port  and  the 
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appropriate  maintenance  base. 

2)  BASEFAILEDAC  will  demand  an  engine  from 
the  port  immediately  in  the  "Continuous 
Pipeline"  configuration  or  at  the  time 
of  its  next  deployment  in  the  "Twisted 
Pipeline"  configuration. 

• If  BASEFAILEDAC  cannot  repair  the  failed 

engine,  l.e.,  BASEFAILEDAC  DESTINBASE,  then: 

1)  Schedule  the  availability  of  the  engine 

at  DESTINBASE  by  invoking  AVAILQ 
procedure  (see  (vi)  below)  with  para- 
meters: (DESTINBASE,  BASEFAILEDAC, 

SHIPNUM (BASEFAILEDAC ) , REPAIRTIME(FAIL- 
DEGREE) ) . 

(Note:  In  the  case  of  a port  which  serves 

as  an  intermediary  point  for  its  carriers, 
provisions  are  made  to  include  the  time 
remaining  in  a carrier's  deployment  by 
using  DUMMYTIME.) 

2)  If  STOCKLEVEL(DESTINBASE)  > 0,  then 
schedule  the  availability  of  a service- 
able engine  of  DESTINBASE  stock  at 
BASEFAILEDAC  by  invoking  AVAILQ  procedure 
(see  (vi)  below)  with  parameters: 
(BASEFAILEDAC,  DESTINBASE,  0,0). 


If  on-hand  stock  does  not  exist  at 
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DESTINBASE,  then  record  a backorder  to 
BASEFAILEDAC  from  DESTINBASE  by  utilizing 
BASEBO  queue.  Queue  backorder  accordingly 
by  invoking  BACKORDQ  procedure  with 
parameter:  (DESTINBASE)  (see  (ix)  below). 

• If  BASEFAILEDAC  can  repair  the  failed  engine, 

i.e.,  BASEFAILEDAC  = DESTINBASE,  then: 
schedule  the  availability  of  this  engine  at 
BASEFAILEDAC  by  Invoking  AVAILQ  procedure  with 
parameters:  (BASEFAILEDAC,  BASEFAILEDAC,  0, 

REPAIRTIME  (FAILDEGREE) ) . 

(li)  ENGAVAILEVENT 

This  procedure  is  invoked  at  the  time  an  engine  in  a 
particular  pipeline  becomes  available  as  on-hand 
serviceable  stock. 

Nomenclature : 

• BASEAVAIL — the  base  number  of  the  base  at  which  an 
engine  has  become  available. 

Synopsis : 

— Determine  BASEAVAIL  from  BASENGAVAIL  queue  and 
rearrange  pointer  mechanism. 

— If  BASEAVAIL  has  a backorder  to  satisfy,  then: 

(The  following  is  performed  utilizing  backorder 
queue  mechanisms.) 

• If  the  backorder  exists  to  BASEAVAIL,  then: 

If  there  are  any  down  aircraft,  then: 


1)  Place  engine  on  first  aircraft  in  down 
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aircraft  queue. 

2)  Determine  the  life  of  the  engine  (flying 
hours ) . 

3)  Determine  the  next  failure  time  of  the 
aircraft . 

If  there  are  no  down  aircraft  at  BASEAVAIL, 
then  place  the  available  engine  in  BASEAVAIL' s 
on-hand  stock. 

• If  the  backorder  exists  to  an  outside  base, 
then  schedule  the  availability  of  the  engine 
at  the  backordered  base  by  Invoklrg  AVAILQ 
procedure  (see  (vi)  below)  with  parameters: 
(backordered  base,  BASEAVAIL,  0,  0). 

— If  there  are  no  backorders  to  satisfy,  place 
available  engine  in  BASEAVAIL's  on-hand  stock. 

(iii)  SHIPEVENT 

This  procedure  is  Invoked  at  the  time  the  state  of  a 
particular  carrier  changes. 

Nomenclature : 

• BASENUM  (I) — base  number  of  carrier  (ship)  number  I_. 

• PORTSCHED  (I,J),  SEASCHED  (I, J)— dock  time, 
deployment  time  of  carrier  (ship)  number  ^ during 

schedule:  J = 1,  2,  • • • , schedule  length. 

(Note:  schedule  may  be  either  deterministic  or 

probabilistic — user  option:  see  Appendix  B. ) 

• SHIPTIME  (I)--tlme  of  next  state  change  of  carrier 
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(ship)  number  I. 

• STATE  (I) — current  state  of  carrier  number 
Synopsis : 

— Determine  which  carrier  is  changing  states  by 
utilizing  SHIPQPT  mechanism. 

— If  this  carrier's  state  change  is  from  'AT  SEA'  to 
'AT  PORT,'  then: 

• Determine  how  long  it  will  be  docked  and  set 
SHIPTIME  accordingly. 

— If  the  carrier’s  state  change  is  from  'AT  PORT'  to 
'AT  SEA, ’ then: 

• Determine  how  long  it  will  be  at  sea  and  set 
SHIPTIME  accordingly. 

• Carrier's  port  attempts  to  restock  carrier 
with  the  number  of  engines  sent  out  during 
last  deployment  plus  any  backorders  from  any 
previous  deployments. 

The  order  of  restocking  is  as  follows: 

1)  Rejuvenate  any  down  aircraft  on  board 
carrier  and  determine  the  life  of  the 
new  engine  plus  the  failure  time  of  the 
aircraft . 

2)  If  there  are  no  down  aircraft,  then 
place  engine  in  carrier's  on-hand  stock. 

If  at  any  time  during  the  restocking  the  port's 
stock  is  depleted,  the  appropriate  number  of 
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backorders  are  recorded  at  the  port, 

— SHIPTIME  (carrier)  is  appropriately  queue  by 
rearranging  SHIPQPT  mechanism  accordingly. 

(Iv)  SAMPLE 

SAMPLE  procedure  is  Invoked  every  SI  simulation  days 
where  SI  is  the  proposed  current  independent  time  length 
(see  Appendix  B) . 

Nomenclature ; 

• NUMSAMPS — Identifies  the  current  sample  number. 

• TOTSAMPLES — total  number  of  samples  desired 
(user  option — see  Appendix  B) . 

• PASSFAIL  (K) — is  1 if  samples  taken  at  iase  K 
passed  the  Independence  test;  is  0 otherwise. 

• PASSBASE — the  number  of  bases  which  passed. 

• TOTNUMACB — the  number  of  bases  which  support  their 
own  flying  activity. 

• PF — the  proportion  of  aircraft  flying  bases  which 
passed  the  Independence  test. 

• SAMPSEATIME  (I,J) — the  amount  of  time  which  carrier 
number  ^ has  spent  in  a deployed  state  during 
sample  number  J. 

Synopsis : 

Its  main  function  is  to  segregate  average  down 
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aircraft  observations  into  Independent  observations.  If 
the  specified  number  of  samples  were  taken,  SAMPLE 
arranges  for  their  Independence  tests. 
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When  SAMPLE  is  called  all  time  distributions  and 
time  related  measurements  are  updated  to  this  point. 

The  sample  number  is  then  advanced. 

If  SAMPLE  is  called  to  terminate  the  final  observa- 
tion CNUMSAMPS  = TOTSAMPLES),  each  base's  down  aircraft 
observations  are  tested  for  independence.  The  results 
of  the  test  are  found  in  the  PASSFALL  array. 

The  proportion  of  passing  bases  is  calculated  and 
then  program  control  is  sent  to  IND_OUT  (see  (xiv)).  If 
control  returns  to  SAMPLE,  more  observations  are  required 
and  further  sampling  is  executed. 

(v)  REPORT 

This  procedure  is  invoked  at  the  termination  of  the 
simulation  run.  Final  statistics  are  reported  to  the 
user. 

(vi)  AVAILQ 

Parameters:  (BASETO,  BASEFROM,  SHIPINDEX,  REPTIME) 

This  procedure  is  Invoked  from  either  PLANEFAILEVENT  or 
ENGAVAILEVENT. 

Synopsis; 

AVAILQ  queues  the  time  at  which  an  engine  becomes 
available  at  BASETO  utilizing  the  ENGAVAILTIME  and 
BASENGAVAIL  queue  mechanisms. 

The  availability  time  may  be  of  the  form: 

• Time  of  engine  availability  = 

Present  time  plus  travel  tire  from  BASEFROM  to 
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BASETO  plus  REPTIME.  REPTIME  can  take  on  any 
positive  value.  For  example,  in  the  case  of 
transporting  a new  engine  from  one  base  to  another, 
REPTIME  = 0.  In  the  case  of  transporting  a failed 
engine  to  be  repaired,  REPTIME  = repair  time  of 
failed  engine. 

• SHIPINDEX  > 0 designates  that  BASEFROM  is  a carrier 
and  is  sending  a failed  engine  to  be  repaired  at 
BASETO.  Since  a carrier  plane  failure  occurs  only 
during  deployment,  an  engine  cannot  be  transported 
to  BASETO  until  the  carrier  has  docked.  Thus, 
time  of  engine  availability  = SHIPTIME  (carrier) 
plus  travel  time  from  the  port  of  the  carrier  to 
BASETO  plus  repair  time  of  the  failed  engine. 

(vii)  FAILGEN 

Parameters;  (PLANEN,  SHIPN,  BASEN) 

This  procedure  is  Invoked  from  PLANEFAILEVENT, 
ENGAVAILEVENT  or  SHIPEVENT. 

Nomenc lature ; 

• DAYSTILFAIL — number  of  days  before  a failed  engine 
causes  PLANEN  to  become  inoperative. 

• FAILTIME— failure  time  of  PLANEN. 

• FLYHRS — flying  hours  of  the  engine  with  the  minimum 
life  on  PLANEN. 

• FOHM  (I,J) — flying  hours  per  month  of  aircraft 
type  I at  base  number  J. 
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• LIFENG  (I,J) — number  of  flying  hours  before  failure 

of  an  engine  In  the  position  on  aircraft  J. 

(I  = 1,  2,  • • • , N:  N E engine  capacity.) 

Synopsis ; 

FAILGEN  Is  called  at  the  time  an  aircraft  Is  removed 
from  the  'down  aircraft*  queue  and  made  operational. 

— Determine  the  position  ^ of  the  engine  with  the 
minimum  life,  l.e.,  FLYHRS  = LIFENG  (1,  PLANEN)  and 
place  this  engine  In  the  first  position. 

— Decrease  the  life  of  all  engines  on  PLANEN  by 
FLYHRS. 

— Determine  DAYSTILFAIL  from  the  expression  FLYHRS/ 
FORM  (PLANEN  type,  BASEN). 

— If  SHIPN  > 0 (BASEN  Is  a carrier),  then  we  have  to 
'pad*  DAYSTILFAIL  so  as  to  make  sure  PLANEN  falls 
during  deployment. 

— Set  PAILTIME  = DAYSTILFAIL  plus  present  time  and 
queue  this  time  accordingly  utilizing  ACFAIL  queue 
mechanism. 

(vlll)  DISTRIBUTIONS 

Parameters:  (BASE,  DIST) 

This  procedure  Is  Invoked  from  either  PLANEFAILEVENT, 
ENGAVAILEVENT,  SHIPEVENT  or  SAMPLE. 

Synopsis : 

DISTRIBUTIONS  Is  called  at  the  time  the  present  stock 
level,  number  of  backorders  or  number  of  down  aircraft 
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is  about  to  change  at  any  base. 

DISTRIBUTIONS  updates  the  weighted  time  frequency 
distribution  and  the  area  under  frequency  distribution 
curve  at  BASE. 

Measurements : 

1)  STOCKDIST  (I,J)  = the  amount  of  time  spent 
at  stock  level  J at  base  I. 


Stock  Level 
at 

Base  I 


Time 


at  time  t,  STOCKDIST  (1,1)  = (t^  - t^)  + 
(ts  - t^)  + (t  - t^^ ) . 

2)  STOCKAREA  (I)  = area  under  the  frequency 
distribution  curve  of  base  1. 

See  above  figure: 
at  time  t,  STOCKAREA  (I) 

= (t^-t^)*!  + + (t2-t2)*l 

+ (tij-t^  ) • 0 + (t-t^^ ) • 1 


A + B + C + D 
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3) 


BODIST  (I,J)  = the 


eimount  of  time  spent  at 


backorder  level  J at  a base 

4)  BOAREA  (I)  = area  under  the  frequency 
distribution  curve  of  base  1. 

5)  DACDIST  (I,J)  = the  amount  of  time  spent  at 
down  aircraft  level  J at  base  ^ which  has 
flying  activity. 

6)  DACAREA  (I,K)  = area  under  the  frequency 
distribution  curve  of  base  1 which  supports 
flying  activity  during  sample  K. 

By  dividing  the  total  simulation  time  in 
equal  samples  and  testing  for  Independence 
among  the  samples,  we  hope  to  perform 
statistical  analysis  on  the  'criterion  value' 
at  each  base  which  is  obtained  directly  from 
DACAREA.  The  test  for  independence  among 
samples  is  crucial  for  obtaining  variances 
and  confidence  intervals  (see  Appendix  D for 
a discussion  on  the  independence  test). 

(Note:  criterion  value  = average  down 

aircraft  at  base  ± / total  number  of  aircraft 
at  base  1.) 
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at  time  DACAREA  (1,1)  = A + B 

DACAREA  (1,2)  = B'  + C + C 
DACAREA  (1,3)  = D 

DIST  is  the  parameter  that  indicates  which  calcula- 
tions are  to  be  performed.  The  values  it  may  take  on 
are  as  follows: 

• •S.TOCK* — indicates  stock  level  measurements  are  to 
be  updated. 

• *B0* — indicates  backorder  measurements  are  to  be 
updated . 

• ’DAC — indicates  down  aircraft  measurements  are  to 
be  updated. 

(ix)  BACKORDQ 

Parameter:  (BASE) 

This  procedure  is  invoked  from  PLANEFAILEVENT  when  a 
backorder  is  recorded  at  BASE. 

BACKORDQ  appropriately  queues  the  backorder. 
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(x)  RNG 


Parameters:  (MEAN,  U) 

This  function  procedure  is  invoked  from  either 
PLANEFAILEVENT,  ENGAVAILEVENT  or  SHIPEVENT. 

Synopsis : 

RNG  returns  an  exponentially  distributed  random 
variable  with  mean  MEAN. 

This  is  achieved  by  using  the  "inverse  transfor- 
mation" technique  of  random  number  generation: 

• Generate  a random  number,  G,  uniformly 
distributed  between  (0,1), 

• Return  the  value  of  F'(U)  where  F*(U)  = Inverse 
of  the  exponential  distribution  function. 


RNG  is  used  to  generate  the  number  of  flying 
hours  until  failure  of  an  engine  on  an  aircraft. 
MTBR  of  the  combination  code  of  the  aircraft,  where 
MTBR  H mean  time  between  removal,  is  the  argument 
which  corresponds  to  MEAN. 

(xl)  SHIPTIME_RNG 

Parameters:  (MU,  SIGMA,  U) 

This  function  procedure  is  invoked  from  SHIPEVENT. 
Synopsis : 

SHIPTIME_RNG  returns  a normally  distributed 
random  variable  with  mean  MU  and  standau’d  deviation 
SIGMA. 
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This  is  achieved  by  using  one  of  the  most 
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efficient  stochastic  variate  generation  techniques 
for  a normally  distributed  random  variable.  This 
technique,  like  most  random  number  generation 
techniques,  is  based  on  the  generation  of  a random 
number,  U,  uniformly  distributed  between  (0,1). 

SHIPTIME_RNG  is  used  when  probabilistic  'dock* 
and  'deployment*  carrier  times  are  desired  (user 
option — see  Appendix  B). 

The  arguments  corresponding  to  MU  and  SIGMA  for 
each  carrier  are  user  Inputs. 

(xli)  SAFETY_STOCK 

Parameters:  (LAMDA,  C) 

This  function  procedure  is  Invoked  once  for  each  base 
during  the  initialization  stage. 

Synopsis ; 

SAFETY_STOCK  returns  the  largest  integer,  S, 
such  that: 

Probability  {x  - S}  - C, 
where  is  defined  as  the  number  of  engines  in 
resupply.  We  assume  that  the  probability  distri- 
bution of  ^ is  Poisson  with  mean  LAMDA.  Thus,  C is 
the  desired  probability  of  having  no  outstanding 
backorders  at  an  arbitrary  point  in  time. 

• LAMDA  = Daily  Demand  Rate  x Average  Resupply 
Time 

• C is  a user  input  for  each  base. 
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SAFETY_STOCK  is  used  In  the  event  a 'DOD- 
requlrement’  calculated  spare  stock  level  is  desired 
instead  of  a user  inputed  level  (user  option — see 
Appendix  B). 

(xiii)  IND_TEST 

Parameters:  (DACAREA,  CRIT_BASE) 

IND_TEST  is  a function  procedure  invoked  from  SAMPLE. 
Nomenclature : 

• XBAR — the  average  of  all  sample  values  taken  at 
the  base. 

• DACAREA — a one  dimensional  vector  in  this  proce- 
dure consisting  of  ordered  sample  observations. 

• CRIT_BASE — the  base  whose  observations  are  being 
tested . 

• SUM  1 — the  statistic  (x^  - 

• SUM  2 — the  statistic  ^ (x^^  - x)^ 

• CN(K) — the  statistic  , SUMl  at  base  K. 

^ ~ 2 ‘SUM 2 


Synopsis : 

IND^TEST  tests  the  observations  at  CRIT_BASE  for 
Independence  and  returns  a value  of  "1"  if  the  series 
is  independent  and  "0"  otherwise. 

IND_TEST  uses  Fishman’s  technique  of  time  series 
independence  testing  [3]  (see  Appendix  D for 
discussion ) . 

(xiv)  IND_OUT 

IND_OUT  is  Invoked  from  SAMPLE,  and  determines 
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whether  enough  bases  have  passed  Independence  tests  (user 
option) 

Nomenclature ; 

• NUMPAILS — the  number  of  times  the  independence 
test  has  failed. 

• LIMTEST — the  user's  limit  on  number  of  tests 
desired . 

Synopsis ; 

If  user  specifications  were  met,  test  results 
are  printed,  listing  number  of  samples,  sample 
interval,  bases,  their  pass/fail  status,  and  its 
test  statistic  calculated  by  IND_TEST.  The  report 
of  entire  simulation  statistics  then  follows. 

Otherwise,  if  more  tests  are  allowed,  successive 
observations  of  downed  aircraft  are  paired  and 
regrouped,  sample  intervals  are  doubled,  and  j 

TOTSAMPLES/2  more  samples  are  regrouped.  Again  test 
results  are  printed  as  in  the  above  paragraph. 

If  the  limit  on  tests  was  reached,  test  results 
are  printed  in  addition  to  some  warnings  and 
suggestions  for  achieving  acceptance  of  Independence. 

A report  is  then  printed. 

(xv)  RES_OUT 

RES_OUT  is  Invoked  from  IND_OUT.  It  simply  prints 
the  results  of  the  independence  test  each  time  the  test 
is  performed.  It  Includes  the  base  name,  a pass/fail 
status,  and  the  CN  test  statistic  for  each  base. 
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Thi’ough  the  specification  of  the  input  values,  many 
options  in  the  system  structure,  performance  and  control  are 
available . 

Input  information  is  in  the  form  of  a stream  of  data 
cards.  Data  can  appear  in  any  column  but  a blank  must  follow 
each  piece  of  data.  Quotation  marxs  are  required  for  some 
data  inputs  (shown  on  example  card).  The  data  cards  must  be 
in  the  order  given  below. 

Instructions  for  coding  along  with  an  example  ’coded’ 
card  are  provided  below: 

Data  Set  #1'.  System  Attributes 
(one  card) 


J52P’ 

’NARP  JAX' 

16  3 

296 

6 

2 

2 

30 

30 

• » 

1 

(a) 

1 

(b) 

1 1 
(c)  (d) 

1 

(e) 

1 

(f) 

1 

(6) 

1 

(h) 

1 

(1) 

1 

(J) 

• k 

(a)  Engine  model  type  for  which  system  is  being  simulated 
(no  more  than  20  characters  long). 

(b)  Depot  name  (no  more  than  20  characters  long). 

(c)  Total  number  of  bases. 

(d)  Total  number  of  carriers,  if  Twisted  Pipeline  Configur- 
ation is  desired: 

"0"  if  Continuous  Pipeline  Configuration  is  desired. 

B1 


B 


B2 


(e)  Total  number  of  aircraft  (number  of  aircraft  at  some 
bases  may  not  be  an  integer  and  should  be  rounded  up 
when  figuring  total  number  of  aircraft  in  system). 

(f)  Total  number  of  aircraft  type, 

(g)  Total  number  of  NAVY  combination  codes  associated  with 
system. 

(h)  The  maximum  engine  capacity  of  all  aircraft  in  system. 

(I)  The  maximum  schedule  length  of  all  carriers  (carrier 
schedules  are  discussed  in  the  explanation  of  Data  Set 
#6),  in  Twisted  Pipeline  Configuration. 

”0”  in  Continuous  Pipeline  Configuration 

(J ) Total  number  of  samples.  This  quantity  is  used  for 
statistical  analysis  purposes  and  must  be  a number 
divisible  by  2.  The  larger  the  value  entered,  the  higher 
the  confidence  in  the  simulation  output  data  along  with 

a higher  cost  associated  with  a longer  computer  run. 

This  quantity  should  be  at  least  20  and  a value  of  30  is 
usually  adequate. 


Data  Set  #2:  Base  Random  Number  Seeds 

(As  many  cards  as  needed) 


09225305  ,165307001 


.059315127  • 


(a) 


(b) 


(k) 


(a) 

Random 

number 

seed 

for 

base 

#1. 

(b) 

Random 

number 

seed 

for 

base 

»2. 

(k)  Random  number  seed  for  base  #k. 


B3 


"0"  and 

k = 

Each  value 
"1"  with  an 

1 2 ^ • 
entered  must  be 
odd  last  digit 

, total  number  of  bases 
a 9-dlglt  number  between 

• 

Data  Set 

#3:  Analysis 

(one  card) 

Parameters 

/ 90 

1 

2 0 
1 1 

1.65  .05 

1 1 

2.045  .05 

1 1 

1 

(a) 

1 1 
(b)  (c) 

1 1 

(d)  (e) 

1 1 
(f)  (g) 

Before  any  statistical  analysis  can  be  performed  on 
the  observed  values  for  each  base's  'criterion*  (=  average 
down  aircraft/ Lotal  number  of  aircraft),  the  variances  of 
these  observed  quantities  m.ust  be  estimated.  The  problem 
arises  from  the  fact  that  the  data  are  serially  correlated 
(that  is,  time  dependent).  In  other  vjords,  if  we  let 
x^(t)  = number  of  down  aircraft  at  base  ^ at  time  then 
x^(t  + s)  could  be  dependent  on  x^(t)  and  x^(t)  could  be 
dependent  on  x^(t  - s)  and  so  on.  The  simple  formulas  for 
computing  sample  variances  are  based  on  independent 
observations.  It  would  seem  reasonable  to  say  that  there 
may  exist  a value  's*  (sometimes  very  large)  such  that 
x^(t  + s)  is  independent  of  x^(t).  If  we  could  find  such  a 
value  's'  and  then  took  observations  *s*  time  units  apart, 
the  values  obtained  would  be  independent  and  the  simple 


1 

] 


i 

i 

i 


1 


variance  formulas  would  be  applicable. 

For  a given  value  of  's,*  SPAERS  tests  the  null 
hypothesis  for  each  base  i,  that  x^(tj^),  x^(t2),  x^Ct^), 

• • • , x^(t^)  is  a random  stream  (l.e.,  Xj^(tj)*s  are 
independent)  where  t^^^  - tj  = s for  J = l,  2,  •••,  N-1 
and  N = total  number  of  samples  specified  on  Data  Set  #1. 

The  test  of  independence  is  said  to  have  'failed’  if  the 
null  hypothesis  is  rejected. 

The  technique  employed  in  testing  the  null  hypothesis 
and  choosing  a new  ’s'  if  the  test  falls  is  outlined  in 
Appendix  D. 


(a)  Initial  length  of  sample  Interval.  This  quantity 
represents  the  user's  guess  to  the  value  of  's.'  If  the 
test  of  independence  using  this  value  of  's’  falls,  then 
's'  is  increased  automatically  in  the  program,  more  data 
are  collected,  and  the  new  's'  is  tested  for  Independence. 
This  procedure  is  repeated  until  the  test  of  independence 
passes  or  the  limit  of  tests  performed  is  reached  (item 
(b)  on  this  card).  An  initial  guess  would  be  to  set  's' 
equal  to  twice  the  maximum  depot  cycle  time  of  all 
bases.  Since  depot  cycle  time  is  typically  the  longest, 
the  system  performance  is  probably  most  affected  by 
depot  activity.  Thus,  by  letting  s = 2 • maximum  depot 
cycle  time,  we  are  allowing  for  the  system  to  'flush' 
Itself  out. 

(b)  The  maximum  number  of  independence  tests  to  be  performed. 
SPAERS  is  constructed  so  that  the  total  number  of  samples 
specified  on  Data  Set  #1  is  constant  no  matter  how  many 
tests  of  independence  are  performed.  For  each  successive 
test  performed,  sample  interval  is  actually  doubled  (see 
Appendix  D).  Thus,  actual  simulation  time  = Run-in 
period  + total  number  of  samples  x Initial  sample  interval 
(item  (a)  above)  x 2^-^  days,  where  k is  the  number  of 
tests  performed.  The  purpose  of  including  this  quantity 
is  to  provide  an  upper  bound  and  hence  a stopping 
criterion  (in  the  Interest  of  computer  time)  if  inde- 
pendence may  not  be  found  within  a reasonable  number  of 
times . 

If  twice  the  depot  cycle  time  is  used  for  the 


initial  sample  Interval,  then  more  than  3 independence 
tests  should  not  be  executed  without  further  investi- 
gation into  the  reasons  for  'high  dependence.' 

In  geneial,  if  a small  initial  sample  interval  is 
chosen  (<  2 x depot  cycle  time),  there  is  a higher  probability 
of  rejecting  the  null  hypothesis  initially,  so  more  inde- 
pendence tests  should  be  allowed. 


(c)  0,  if  all  bases  are  required  to  pass  independence  test. 

Base  number  of  a 'critical'  base  which  is  believed  to  be 
the  highest  candidate  for  failure  of  independence.  In 
choosing  a 'critical'  base  it  is  the  belief  that  if  this 
base  passes  the  test  for  Independence,  then  it  is 
reasonable  that  the  remaining  bases  would  also  pass  the 
test.  Thus,  it  is  only  necessary  to  require  that  the 
independence  test  pass  at  the  'critical'  base. 

(d)  and  (e)  Parameters  for  the  test  of  independence. 

Let  ^ denote  the  value  entered  in  (d)  and  a the  value 
entered  in  (e ) . 

Then,  ^ is  the  1 - a significance  point  of  the  standard 
normal  distribution. 

Graphically : 


^ is  the  test  statistic  and  a is  the  significance 
level  for  testing  the  iiull  hypothesis,  that  x(t,), 

x(t2),  * • * > x^t^)  are  Independent. 

‘^Flrst,  the  user  must  choose  a significance  level  a. 
We  note  that: 

a = Probability  (rejecting  Hq  when  Hq  Is  true).  A 
small  a level  means  that  there  is  a higher  probability 
of  accepting  Hq  when  it  is  false  and  would  lead  to  an 
underestimation  of  the  variance  of  On  the  other 

hand,  a large  a level  means  that  there  is  a higher 
probability  of  rejecting  Hq  when  it  is  true  and  would 
lead  to  a less  efficient  estimate  of  the  confidence 


Interval  of  Significance  levels  of  0.05  and  0.10 

are  typically  used.  Levels  smaller  than  0.025  and 
larger  than  0.25  should  be  avoided. 

Z is  obtained  from  standard  normal  tables  such  that; 
Probability  (standard  normal  - Z)  = 1 - a. 

It  would  seem  reasonable  to  accept  the  hypothesis  of 
independence  for  all  bases  if;  the  number  of  bases  that 
failed  the  test  of  Independence  divided  by  the  total 
number  of  bases  tested  is  less  than  a,  since  there  is  an 
error  *a'  of  rejecting  Hq.  This  is  in  fact  the  passing 
criterion  for  system  Independence  used  by  SPAERS. 

(f)  and  (g)  Confidence  Interval  parameters. 

Let  ^ denote  the  value  entered  in  (f)  and  y the  value 
entered  in  (g). 

^ and  Y ai’e  used  to  compute  confidence  intervals  on  the 
mean  criterion  at  each  base. 

Interpretation  of  ^ and  y is  as  follows: 

(1)  Pr(L^  ^ U ^ L2>  = 1 - Y, 

where:  L,  = x - — 

^ /n 

L.  = 3T  + ^ 

W = true  mean 

and,  X = sample  (observed  mean) 

s = sample  standard  deviation 
n = number  of  samples  observed 

Y 

t = 100(1  - percentile  of 
the  student's  t-dlstribu- 
tlon  with  n - 1 degrees  of 
freedom. ** 

In  short,  equation  (1)  defines  limits  L,  and  L2  such 
that  there  is  100  (1  - y)%  confidence  that  the  true  mean 
lies  between  these  bounds. 

Commonly  used  values  for  y are  0,05  and  0.01  to  give 


• Reference  [3],  page  7. 

**  Referei.ee  [4], 


r 
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a "confidence  Interval"  of  95!^  and  99/S,  respectively. 

Once  a y level  has  been  chosen,  a £ value  may  be 
obtained  from  t-distribution  tables.  t is  chosen  such 
that ; 


Pr{x^t}=l-^ 


where  x is  a student's  t-distributed  random  variable  with 
n - 1 degrees  of  freedom.  (n  = total  number  samples; 
entered  on  Data  Set  #1.) 


Data  Set  #4:  System  Combination  Codes 

(As  many  cards  as  needed) 


102 

1 

104 

1 

• • • 

108 

1 

• • • 

1 

(a) 

1 

(b) 

• • • 

1 

Ck) 

• • • 

(a)  Navy  combination  code.  Throughout  the  remainder  of  this 
section  and  in  the  simulation,  CC102  will  be  indexed  by 
a "1." 

(b)  Navy  combination  code.  Throughout  the  remainder  of  this 
section  and  in  the  simulation,  CC104  will  be  indexed  by 
a "2." 


(k)  Navy  combination  code.  Throughout  the  remainder  of  this 
section  and  in  the  simulation,  CC108  will  be  indexed  by 
a "k." 


k - 1,  2,  3, 


• • , total  number  of  combination  codes 


in  system 


f 


Data  Set  #5’.  Aircraft  Type  Attributes 
(one  card  for  each  aircraft  type) 


/ *AGE»  2 2 

I I I 

(a)  (b)  (c) 

(a)  Navy  aircraft  type  name  (no  more  than  7 characters  long). 

(b)  Engine  capacity. 

(c)  Combination  Code  index  (see  Data  Set  #k) , In  this  case 
the  combination  code  of  AGE  is  104. 

For  the  remainder  of  this  section  and  in  the  simu- 
lation, an  aircraft  type  will  be  Indexed  by  the  order  number 
in  which  it  was  ’inputed.’ 


Data  Set  #S:  Carrier  Attributes 

(one  card(s)  for  each  carrier) 


15  20  60  10  :50  10  15 


(f) 


(a)  Mean  'docking'  time  (time  spent  in  port)  in  days,  if 
probabilistic  times  are  desired; 

0,  if  deterministic  'docking'  times  are  desired. 

(b)  Mean  'deployment'  time  (time  spent  at  sea)  in  days,  if 
probabilistic  time  are  desired; 

0,  if  deterministic  'deployment'  times  are  desired. 


(c)  standard  deviation  (in  days)  of  docking  and  deployment 
time'  is  probabilistic  times  are  desired; 

0,  if  deterministic  docking  and  deployment  times  are 
desired . 

(d)  0,  if  probabilistic  carrier  times  are  desired.  Docking 
and  deployment  times  will  thus  be  random  variables, 
normally  distributed  with  mean  (a)  and  (b)  and  standard 
deviation  (c). 

1,  if  deterministic  schedules  are  desired. 

(e)  30,  if  (d)  is  equal  to  0; 

The  desired  number  of  times  before  deterministic  schedule 
is  repeated,  if  (d)  is  equal  to  1. 

The  two  numbers  in  (e)  have  different  meaning  and 
should  not  be  confused.  In  the  probabilistic  case,  we 
need  to  know  the  future  docking  and  deployment  times  when 
determining  a failure  time  for  an  aircraft  since  a 
failure  during  the  ’docking’  state  is  not  permitted 
(aircraft  does  not  fly).  Thus,  the  number  entered  in 
this  case  Indicates  that  we  may  be  safe  in  not  scheduling 
a failure  during  the  ’docking’  state  if  this  number  of 
future  docking  and  deployment  times  are  known.  The 
number  30  is  sufficient  for  MTBR  (mean  time  between 
removal)  = 226  but  should  be  increased  if  a MTBR  much 
larger  than  226  is  used. 

On  the  other  hand,  the  deterministic  schedule 
repeats  in  cycles  and  thus  knowledge  of  future  carrier 
times  is  present  at  all  times.  The  number  entered  here 
thus,  designates  the  schedule  cycle  length  or  in  other 
words,  the  number  of  schedules  in  the  deployment  (docking) 
cycle . 

A probabilistic  schedule  may  be  desired  when  a 
particular  carrier  has  unusual  but  consistent  and/or 
required  deployment  and  docking  times,  as  may  be  the 
case  during  wartime. 

(f)  a blank,  i.e.,  no  numbers  entered,  if  (d)  equals  0; 

Docking  time  of  first  schedule. 

Deployment  time  of  first  schedule. 


Docking  time  of  ^th  schedule. 

Deployment  time  if  ^th  schedule,  v/here  i = 1,  2,  • • • , 
schedule  length  (value  entered  in  (e)),  if  (d)  equals  1. 
(Note;  docking  time  = time  spent  in  port;  deployment 
time  = time  spent  at  sea.) 


BIO 

The  order  of  the  cards  of  Data  Set  ^6  should  be  the 
same  as  the  ship  numbers  of  the  carriers  (see  Data  Set  #10 
for  ship  numbers). 

(Important ; Data  Set  #6  Is  omitted  If  Item  (d)  on  Data 
Set  #1  equals  0,  l.e..  Continuous  Pipeline 
Configuration. ) 

Data  Set  #7:  Mean  Time  Between  Removal  (MTBR’s) 

(As  many  cards  as  needed) 


187  226  ...  3212 


(a)  (b)  ...  (k) 


(a)  MTBR  as::-oclated  with  the  combination  code  Indexed  by  1 
(see  Card  Set  #4). 

(b)  MTBR  associated  with  the  combination  code  Indexed  by  2. 
(k)  MTBR  associated  with  the  combination  code  Indexed  by  k. 


k = 1.  2 


total  number  of  combination  codes 
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Data  Set  #8:  Combination  Code  Maintenance 

Probabilities 

(one  card  for  each  combination  code) 


.116  .362  .230  .115  .117 

I I I I I 

(a)  (b)  (c)  (d)  (e) 

Probability  that,  given  an  engine  failure,  the 
degree  of  maintenance  is: 

(a)  3rd  degree 

(b)  2nd  degree 

(c)  1st  degree 

(d)  (-l)st  degree  (maintenance  performed  at  depot  on  a 1st 
degree  failed  engine) 

(e)  overhaul  or  0^^  degree. 

The  order  of  the  cards  of  Data  Set  #8  should 
correspond  to  the  index  of  the  respective  combination  code(s)o 


Data  Set  Engine  Repair  Times 

(one  card) 


T 

19 

1 

2iJ 

1 

23 

1 

30 

1 

(a) 

1 

(b) 

1 

(c) 

1 

(d) 

1 

(e) 

The  repair  time  (in  days)  required  for  the  following 


.i 


. i 


! 

I 


i 


\ 


} 

{ 


r 
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degree  maintenance: 

(a)  3rd  degree 

(b)  2nd  degree 

(c)  1st  degree 

(d)  (-l)st  degree 

(e)  overhaul  or  0^^  degree. 


Data  Set  ^10:  Base  Synopsis 

(one  data  set  for  each  base) 

Data  Set  #10A:  Base  Attributes 

(one  card) 


'CHERRY  POINT'  40020131. 01  2 


(a) 


(b)  (c)  (d)  (e)  (f)  (g)  (h)  (i)  (j)  (k) 


(a)  Base  name  (cannot  exceed  20  characters  in  length). 


(b)  Base  number  or  the  order  the  base's  data  set.  In  this 
example,  CHERRY  POINT  is  the  4th  base  to  be  read. 

(c)  Carrier  number  or  the  order  in  which  it  is  read  in 
comparison  to  other  carriers,  if  the  base  is  a carrier 
and  Twisted  Pipeline  Configuration  is  desired; 

0,  if  the  base  is  not  a carrier  or  if  the  base  is  a 
carrier  but  a Continuous  Pipeline  Configuration  is 
desired . 

(d)  1,  if  base  is  an  'acting'  port; 

0, 


otherwise 
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'Acting'  port  is  defined  as  a port  which  serves  as 
an  intermediary  point  for  all  engine  transactions  between 
its  carriers  and  the  appropriate  higher  echelon  maintenance 
bases . 

This  concept  was  Introduced  to  allow  for  the 
inclusion  of  a base  which  neither  supports  flying  activity 
nor  maintenance  activity  but  serves  as  a port  to  one  or 
several  carriers.  This  type  of  base,  in  many  cases,  is 
currently  ignored  in  the  analytical  models.  The  reason  for 
this  is  that  the  pipeline  length  from  a carrier  to  base  i, 
say,  is  equal  to  the  sum  of  the  pipeline  length  from  a carrier 
to  its  port  and  the  pipeline  length  from  that  port  to  base 
The  pipeline  lengths  in  a continuous  pipeline  configuration 
remain  constant  and  since  this  port  supports  neither  flying 
nor  maintenance  activity,  it  serves  no  computational  purpose. 
On  the  other  hand,  in  the  real  world  system  this  type  of 
port  does  serve  as  an  intermediary  point  and  although  its 
exclusion  may  have  no  apparent  effect  in  the  analytical 
model,  the  potential  hazards  in  a twisted  pipeline  configur- 
ation are  worth  noting. 

Excluding  this  port  in  a twisted  pipeline  configur- 
ation would  mean  to  consider  the  carrier's  next  echelon 
maintenance  base  as  its  port.  This  would  mean  that  this 
'artificial'  base  has  accessibility  to  the  carrier  engines 
that  were  in  repair  at  the  time  of  a carrier  departure  from 
port  but  have  become  available  as  serviceable  on-hand  stock 
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since.  This  accessibility  would  last  until  the  time  of  the 
carrier's  next  deployment  (i.e.,  the  time  remaining  in 

present  deployment  plus  the  time  spent  in  port).  Further- 
more, the  fact  that  this  base  employs  the  carrier  engines 
increases  the  probability  that  the  carrier  will  depart  with 
backorders.  These  factors  will,  of  course,  affect  the 
performance  of  the  system,  especially  if  this  base  serves 
as  the  maintenance  activity  for  many  other  bases. 

Even  in  the  continuous  pipeline  configuration,  the 
artificial  port  should  be  avoided,  since  increasing  the  spare 
stock  level  at  this  port  to  compensate  for  poor  carrier 
performance  would  have  similar  effects  to  that  of  the  twisted 
pipeline  configuration:  this  port  would  have  accessibility 

to  these  engines  when  not  in  demand  from  the  carrier. 

It  is  the  belief  of  the  authors  that  all  "real" 
ports  should  be  included  whether  continuous  or  twisted 
pipeline  configurations  are  desired  since  their  inclusion  is 
more  representative  of  the  real-word  system.  These  "real" 
ports  would,  as  noted,  have  a "1"  under  item  (d)  (a  "0" 
would  be  entered  for  artificial  ports). 

(e)  number  of  aircraft  type. 

(f)  base  number  of  the  base  which  serves  as  a port  for  this 
base,  if  this  base  is  a carrier  and  Twisted  Pipeline 
Configuration  is  being  employed; 

0,  if  base  is  not  a carrier  or  Continuous  Pipeline 
Configuration  is  being  employed. 


w 


(g)  0,  if  it  is  desired  to  have  the  prograjn  compute  spare 
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stock  level  for  this  base,  using  the  DODI  ^230.4 
calculation; 

1,  if  it  is  desired  to  'input*  the  spare  stock  level. 

(h)  0,  if  (g),  above,  is  at  level  0. 

Desired  spare  stock  level,  if  (g),  above,  is  at  level  1. 

(i)  A desired  probability  of  having  no  outstanding  backorders 
at  an  arbitrary  point  in  time,  if  (g),  above,  is  at 
level  0. 

0,  if  (g),  above,  is  at  level  1. 

(J  ) Desired  criterion. 

(k)  Weight  of  importance  relative  to  other  bases.  (A  desired 
criterion  of  0.01  at  base  jL  and  0.02  at  base  would  mean 
that  base  is  seen  as  being  twice  as  important  as  base 
and  consequently  the  weight  of  base  would  be  twice  that 
of  base  j_.  ) 

Data  Set  #10B:  Base  Maintenance  Activity 

(one  card) 


7 7 6 8 8 

I I I I I 

(a)  (b)  (c)  (d)  (e) 

Maintenance  activity  (base  number)  of  this  base's 
engine  failures  of  degree: 

(a)  3 

(b)  2 

(c)  1 

(d)  (-1)  (always  the  same  as  (e)) 

(e)  overhaul  or  0 


( Important ; 1)  For  a carrier  which  is  associated  with  a 

base  having  ’acting*  port  indicator  equal 
to  1 (see  Data  Set  if'lOA,  (d)),  the  values 
entered  for  degree  failures  unable  to  be 
maintained  by  the  carrier  should  cor’^’espond 
to  the  base  number  of  its  Acting’  port  (l.e., 
the  base  with  ’acting’  port  indicator  equal 
to  1).  Similarly,  the  values  entered  for 
degree  failures  of  a base  with  ’acting’  port 
Indicator  equal  to  1 should  correspond  to  the 
base  number  of  the  appropriate  carrier 
maintenance  base. 


2)  For  a base  which  does  not  support  flying 
activity,  its  base  nuraber  should  be  entered 
under  its  highest  repair  degree  and  "0" 
everywhere  else.) 


Data  Set  #10C:  Base  Travel  Times 

(As  many  cards  as  needed) 


0 0 0 10  3 


(a)  (b)  (c)  (d)  (e) 


Travel  time  (in  days)  to: 
(a)  Base  #1 


B17 


(b)  Base  §2 

(c)  Base  #3 

(d)  Base 

(e)  Base  #5 

(k)  Base  #k 


K = 1,  2,  3,  * • • , total  number  of  bases. 

(Note;  Travel  time  from  base  ^ to  base  ^ is, 
of  course,  equal  to  0. 

All  carrier  in  Twisted  Pipeline  Configuration 
have  0 travel  time  to  all  bases  in  system.) 


Data  Set  #10D:  Base  Aircraft  Type  Attributes 

(one  card  for  each  aircraft  type) 


27.29 

1 

(b) 


(c) 


(a)  The  order  number  or  index  of  aircraft  type  (see  Data 
Set  #5). 

(b)  Flying  hours  per  month  per  aircraft  of  this  type. 

(Note:  Given  the  fact  that  carriers  in  Twisted  Pipeline 

Configuration  may  only  support  flying  activity 
during  deployment,  this  value  is  adjusted  in  the 
program  such  that  the  average  total  flying 


B18 


activity  in  both  the  Twisted  and  Continuous 
Pipeline  Configurations  are  the  same.) 

(c)  Number  of  aircraft  of  this  type  at  base. 

(Non-integer  values  should  not  be  rounded  up.  The 
program  performs  the  'rounding  up'  of  number  of  aircraft 
and  adjusts  the  flying  activity  such  that  total  flying 
activity  is  the  same.) 


Summary 


Data  Set  #1  - 
Data  Set  #2  - 

Data  Set  #3  - 
Data  Set  - 

Data  Set  #5  - 

Data  Set  ^6  - 

Data  Set  #1  - 

Data  Set  ^8  - 

Data  Set  H9  - 
Data  Set  #10- 


System  Attributes  (one  card) 

Random  Number  Seeds  (one  for  each  base;  as 
many  cards  as  needed) 

Analysis  Parameters  (one  card) 

System  Combination  Codes  (as  many  cards  as 
needed) 

Aircraft  Type  Attributes  (one  card  for  each 
aircraft  type  in  system) 

Carrier  Attributes  (one  card  for  each  carrier; 
this  set  is  omitted  in  Continuous  Pipeline 
Configurations ) 

MTBR's  (one  for  each  combination  code--ln 
same  order  as  Data  Set  #4;  as  many  cards  as 
needed ) 

Combination  Code  Maintenance  Probabilities 
(one  card  for  each  combination  code — same 
order  as  Data  Set  §k) 

Engine  Repair  Times  (one  card) 

Base  Synopsis  (one  data  set  for  each  base) 


r 
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••  Data  Set  i^lOA  - Base  Attributes  (one  card) 

••  Data  Set  iS<10B  - Base  Maintenance  Activity  (one 

card ) 

••  Data  Set  #10C  - Base  Travel  Times  (as  many  cards 

as  needed) 

••  Data  Set  )5*10D  - Base  Aircraft  Type  Attributes 

(one  card  for  each  aircraft  type; 
data  set  omitted  if  base  does  not 
support  flying  activity) 


Appendix  C;  SPAERS  Output 


Contents 


Cl.  Discussion  of  Performance 

Measures  

CII.  Example  Listing  of  SPAERS  Output  . . 


Cl 
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Cl.  Performance  Measures:  Data  Output 


At  the  end  of  the  simulation  run  statistics  summar- 
izing the  performance  of  each  base  and  the  overall  system 
are  reported  to  the  user. 

A discussion  of  these  performance  measures  follows: 

* Rate  is  the  probability  that  no  backorders  exist  at 

any  one  point  in  time.  If  we  let  NI  represent  net 
inventory:  that  is,  on-hand  stock — backorders,  then 

ready  rate  is  given  by 

Probability  {NI  ^ 0}. 

Expressed  in  another  way,  ready  rate  is  the  proportion 
of  time  the  base  has  spent  in  a state  of  *0  backorders,’ 
or  equivalently,  the  proportion  of  time  spent  in  a state 
of  ’net  inventory  - 0.’ 

* Rate  is  the  proportion  of  demands  that  can  be  satis- 
fied immediately.  Equivalently,  fill  rate  is  the 
probability  that  an  engine  is  available  when  a demand 
occurs.  Thus,  fill  rate  is  given  by 

Probability  {NI  > 0}. 

We  note  that  Ready  Rate  = Fill  Rate  + Pr{NI  = 0}. 
Expressed  in  another  way,  fill  rate  is  the  proportion  of 
time  the  base  has  spent  in  a state  of  ’net  inventory  - 1. 
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• Immediate  Satisfaction  Rate  Is  the  proportion  of  demands 

that  were  satisfied  Immediately  during  the  simulation 
run.  Immediate  satisfaction  rate  Is  a measure  of 
frequency:  total  number  of  demands  satisfied 

immedlately/total  number  of  demands.  Fill  Rate  Is  a 
measure  of  time. 

Note  that  as  time  -*■  Immediate  Satisfaction 
Rate  -*■  Fill  Rate. 

• Average  On- Hand  Stock 

This  quantity  Is  obtained  from  .(^ase  1_) 

^ ’’  simulation  time 

(see  section  B.  (2)(vlll)). 

• Average  Backorders 

Similarly,  this  quantity  Is  obtained  from 

BOAREA  (base  1) 
simulation  time  * 

• Average  Backorder  Time  is  defined  as  the  average  time  until 

a backorder  Is  satisfied. 

• Average  Down  Aircraft 

This  quantity  is  obtained  from  A I 

^ ^ simulation  time 

• Average  Down  Aircraft  Time  is  defined  as  the  average  time 

a down  aircraft  remained  Inoperative, 

j average  down  aircraft 

• Criterion  is  defined  as  — • 

total  number  of  aircraft 

The  standard  deviation  along  with  a (1  - <x)%  confidence 

Interval  are  reported. 

(a  Is  a user  Input — see  Appendix  B. ) 

• Non-Zero  Distribution  Values  for  net  Inventory  and  down 


aircraft  are  reported.  Distribution  value  of  a variable 


is  defined  here  as  the  proportion  of  time  the  variable 
spent  at  some  level 

For  a base  which  repairs  only  its  own  engines,  the  back- 
order and  down  aircraft  statistics  are  equivalent. 

Down  aircraft  statistics  are  omitted  for  bases  with  no 
'flying  activity.' 

Special  provisions  must  be  made  for  the  reporting  of  a 
base  which  serves  as  a port  in  the  Twisted  Pipeline 
Configuration.  The  two  events:  backorders  > 0 and 

on-hand  stock  > 0,  are  mutually  exclusive  for  all  other 
bases  since  backorders  must  be  satisfied  at  the  time  an 
engine  becomes  available.  In  the  case  of  a port,  the 
carrier  for  which  it  is  backordered  to  may  be  in  deploy- 
ment state  when  an  engine  becomes  available  and  thus, 
the  possibility  of  on-hand  stock  > 0 and  backorders  > 0. 

Because  of  this  characteristic,  it  would  not  make 
sense  to  talk  about  net  inventory  at  a port  since  'net 
Inventory*  = 0 could  mean  on-hand  stock  = 1 and  backorders 
= 1 , or  on-hand  stock  = 2 and  backorders  = 2,  etc.  Thus, 
non-zero  distribution  values  are  reported  for  'on-hand 
stock'  and  'backorders'  separately. 

Furthermore,  the  'ready  rate'  and  'fill  rate'  are 
omitted  since  these  measures  have  no  meaning  under  the 
'unique'  conditions  of  a port  in  the  Tv.’isted  Pipeline 
Configuration. 
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System  Performance 

The  criterion  value  obtained  from  the  simulation 
run  is  compared  to  the  desired  criterion  (user  input 
see  Appendix  B)  for  each  base  with  flying  activity. 
Furthermore,  two  system  measures  are  reported: 

(1)  System  Performance  is  defined  as: 

K 

^ average  down  aircraft  at  base  i 

i=l , 

K 

I number  of  aircraft  at  base  i 
i=l 

where  k = number  of  bases  with  flying 
activity. 

(2)  Weighted  Average  Performance  is  defined  as: 

k 

J W.  • Criterion  of  base  i 
i=l  ^ 

k 


where  = weight  of  importance  of  base  i 
(user  input — see  Appendix  B) . 


CII.  Example  Listing  of  SPAERS  Output 


The  format  of  SPAERS  output  is  as  follows: 

First,  a report  of  system  parameters  is  printed. 

This  serves  as  a check  for  user  input  errors. 

The  spare  stock  level  for  each  base  along  with  an 
indicator  to  whether  this  level  was  an  input  or  a DODI-^230^ 
calculation  is  reported  next. 

A synopsis  of  each  test  for  independence  follows. 

The  synopsis  includes  the  number  of  samples,  the  sample 
Interval,  the  pass/fall  status  for  each  base  and  the  value 
of  its  test  statistic  for  each  test  performed.  If  the  limit 
on  tests  is  reached,  some  warnings  and  suggestions  for 
achieving  acceptance  of  Independence  are  also  printed. 

Finally,  performance  measures  for  each  base  followed 
by  system  measures  are  reported. 

An  example  of  SPAERS  output  follows. 
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The  technique  employed  by  SPAERS  for  deriving 


Interval  estimates  of  sample  means  from  the  output  of  a 
simulation  run  was  presented  by  George  S.  Fishman  in  a 
technical  report  in  August  1975.*  The  development  of  the 
method  is  based  upon  the  concept  of  'batch  means.'  A sample 
sequence  is  divided  into  batches  of  equal  size  from  which  a 
sample  mean  in  each  batch  is  computed  along  with  an  estimate 
of  the  variance  of  the  grand  sample  mean  over  all  batches. 

The  underlying  assumption  is  that  the  batch  means  are 
independent.  The  technique  uses  the  von  Neumann  ration**  to 
test  for  independence  among  batches. 

Let  x^,  x^,  • • • , x^  be  a sequence  of  identically 
distributed  random  variables.  We  want  to  test  the  hypothesis 
that  the  x^  are  independent.  Consider  the  statistic:** 

C = 1 - 

" r n 2-| 

'■‘I  ■ 

-1  r 

where  x„  = n Z 

" i=l  ^ 

* Reference  [3]. 

**  Reference  [7]. 
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If  the  are  independent,  then 

E(C^)  = 0. 

Furthermore,  if  the  data  are  also  from  a normal  population, 
then 

Certain  properties  show  that  asymptotically  has 
the  same  "principal  distributional  characteristics"*  as  in 
the  normal  case  (i.e.,  even  if  the  are  not  from  a normal 
population).  can,  thus,  be  treated  as  though  the  x^  were 

normal,  for  large  n.  Furthermore,  becuase  of  this  good 
approximation  of  the  normal  distribution  to  the  distribution 
of  in  the  normal  case,  C^//  (n-2)/(n^  - 1)  may  be  treated 
as  a standard  normal  (that  is,  normal  with  mean  = 0,  standard 
deviation  =1)  when  n is  large.  For  the  0.05  significance 
level  and  n - 8,  this  approximation  introduces  an  error  of 
less  than  0.2^.** 

Thus,  may  be  regarded  as  being  symmetrically 
distributed  about  0 (in  samples  from  a normal  population). 
Significant  departure  of  from  0 may  be  taken  as  indicating 
non-randomness  in  the  time  series.  In  general,  positive 
values  of  indicate  the  presence  of  positive  correlation 
while  negative  values  Indicate  negative  correlation.*** 


* Reference  [3]. 

**  Reference  [8]. 

***  Reference  [9]- 
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The  test  of  hypothesis  then  is;  Accept  where 

^o'  ^1’  ^2*  ’ * * * Independent,  if  C^//(n-2 )/(n2-l ) - 

Z^,  where  Is  the  1 - a significance  point  of  the  standard 
normal  distribution.  Here  a denotes  the  Type  I error:  a = 

Probability  of  rejecting  v;hen  is  true.  The  reason  why 
the  test  is  one-sided  is  that  only  positive  correlation  is 
typically  present  in  simulation  time  series  of  the  kind 
found  in  SPAERS.  That  is,  if  many  'down'  aircraft  are 
observed  at  time  t,  one  would  not  expect  to  see  few  or  no 
'down'  aircraft  at  time  t+e  and  again  at  time  t+2e  (where  e is 
some  small  time  interval).  Rather,  if  many  'down'  aircraft 
are  observed  at  time  then  many  'down'  aircraft  will 
probably  be  observed  at  time  t + e and  time  t + 2e.  Thus, 
negative  values  of  will  not  be  significant  to  reject  the 
null  hypothesis o 

The  original  study  of  involved  a single  test. 

George  S.  Fishman  outlined  in  his  technical  report  [3]  an 
iterative  procedure.  It  is  this  procedure  which  SPAERS 
employs . 

The  set  up  is  as  follows: 


1 1 

?:"■  j 

Period 

Sample 

1 

Sample 

Z 

• • • 

Sample 

i 

t « • 

Sample 

n 

K- 

T 

T H 

h* T H 

. T - - i 

|<- SIMULATION  TIME  ^ 


F 


DiJ 


^ / y(t)dt  1 = 1,  2,  •••,  n 

(i-l)T 

where  Y(t)  = number  of  down  aircraft  at  time  ^ 

n =»  total  number  of  samples 

= time  weighted  average  down 
aircraft  of  sample  ^ 

The  estimate  of  average  down  aircraft  over  the  total 

simulation  time: 

_ n , nT 

X = n“  I X.  = -ijr  / Y(t)dt 
1=1  o 

If  x^’s  are  Independent,  then  the  variance  estimate 

= [n(n  - 1)]"^  I (x.  - . 

1 = 1 ^ 

To  test  if  the  x^  are  independent,  compute  and 

accept  the  hypothesis  of  Independence  if  (n-2 )/ (n^-i ) 1 Z^. 

If  the  hypothesis  is  rejected,  then  create  a new 

sequence  x' • ^*n/2  averaging  pairs  of 

non-overlapping  successive  samples.  That  is:  Let 

x’l  = (Xj^  + X2)/2,  x*2  = (x^  + x^)/2,  • • • , 

x'^^2  " ^^n-1  have  n/2  samples  of  length 

2T  and  thus  n/2  observations.  Test  the  new  C . If  the 

n 

hypothesis  is  again  rejected,  then  create  a new  sequence 
x"i,  x"^,  • * • , ^"n/4  average  pairs  of  non-overlapping 
successive  samples.  We  now  have  n/^'  observations.  As  the 


procedure  continues,  the  data  are  effectively  batched  In 
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groups  of  1,  2,  4,  8,  l6,  . , , non-overlapping  observations. 
As  the  batches  increase  in  size,  one  would  expect  the 
correlation  between  the  observations  to  decrease.  Eventually 
the  null  hypothesis  is  accepted.  When  it  is  accepted,  the 
sample  variance  is  computed  and  a confidence  Interval  for  3T 
may  be  determined  using  the  t-dlstrlbutlon.  (It  should  be 
noted  that  n in  this  case  must  be  a power  of  2.) 

One  problem  with  this  iterative  approach  is  that 
the  number  of  observations  decreases  rapidly  as  the  number  of 
tests  that  have  to  be  performed  increases.  The  fewer  the 
observations,  the  less  confident  we  are  in  our  observed 
values.  To  alleviate  this  potential  problem,  the  number  of 
samples,  n (observations),  remains  constant  in  SPAERS  no 
matter  how  many  tests  are  performed.  What  is  dependent  on 
the  number  of  tests  performed  is  simulation  time.  That  is, 
if  the  first  test  of  Independence  fails  (reject  null 
hypothesis),  then  create  a new  sequence  x’j*  ^'2*  * * ** 
x’n/2»  before.  The  time  interval  of  x'^^^  is  2T  and  we 
have  n/2  observations.  At  this  point,  the  simulation  is 
continued  and  n/2  more  observations,  2T  time  units  apart,  are 
collected.  Thus,  at  the  time  the  next  test  of  Independence 
is  performed,  we  will  have  n observations.  This  procedure 
continues  until  the  test  of  Independence  passes  or  the  limit 
of  tests  performed  is  reached  (user  input — see  Appendix  B). 
(It  should  be  noted  that  the  only  restriction  on  n,  total 
number  of  samples,  is  that  it  be  divisible  by  2.) 
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